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Abstract 



The upper limit of the total cross section for quasi-free pn — > pnv[ reaction has been 
determined in the excess energy range near the kinematical threshold. 
The measurement has been carried out at the COSY-1 1 detection setup using a pro- 
ton beam and a deuteron cluster target. The identification of the T|' meson has been 
performed using the missing mass technique. The energy dependence of the up- 
per limit of the cross section is extracted with a fixed proton beam momentum of 
Pbeam = 3.35 GeV/c and exploiting the Fermi momenta of nucleons inside the deuteron. 
The data cover a range of centre-of-mass excess energies from to 24 MeV. 
The experimentally determined upper limit of the ratio = cipp^pprf) ' wn i cn is 
smaller than the ratio for the T] meson, excludes the excitation of the Sn(1535) reso- 
nance as a dominant production mechanism of the r\' meson in nucleon-nucleon col- 
lisions. At the same time, the determined upper limits of go in the direction of 
what one would expect in the glue production and production via mesonic currents. 
For quantitative tests of these mechanisms an order of magnitude larger statistics and a 
larger energy range would be required. This can be reached with the WASA-at-COSY 
facility. 

Streszczenie 



Gorna granica calkowitego przekroju czynnego dla kwazi-swobodnej reakcji 
pn — > pnr\' zostala wyznaczona w przyprogowym obszarze energii. 
Pomiar zostal przeprowadzony wykorzystujac uklad detekcyjny COSY-1 1, wiazke 
protonowa oraz klastrowa tarcze deuteronowa. Mezon T|' zostal zidentyfikowany 
przy uzyciu techniki masy brakujacej. Zaleznosc gornej granicy calkowitego 
przekroju czynnego od energii (funkcje wzbudzenia) przy stalej wartosci pedu wiazki 
Pbeam = 3.35 GeV/c, uzyskano dzieki pedowi Fermiego nukleonow w deuteronie. 
Dane zostaly zmierzone w przedziale energii nadprogowej Q od MeV do 24 MeV. 
Otrzymana gorna granica stosunku R^/ calkowitych przekrojow czynnych dla reakcji 
pn — > pnr\' i pp — > ppr\', ktorej wartosc okazala sie bye znacznie mniejsza od analog- 
icznego stosunku zmierzonego dla mezonu T|, wyklucza hipoteze iz dominujacym pro- 
cesem w produkeji mezonu r\' jest wzbudzenie rezonansu Sn(1535). Rownoczesnie, 
uzyskany wynik nie wyklucza innych mechanizmow produkeji mezonu T|' takich jak 
prady mezonowe czy wzbudzenie gluonow. 

Dla przeprowadzenia ilosciowych testow mozliwych mechanizmow produkeji mezonu 
T|' wymagane jest wykonanie eksperymentu z wieksza statystyka oraz w wiekszym za- 
kresie energii. Taka mozliwosc daje uklad detekcyjny WASA-at-COSY. 
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I. Introduction 



The main aim of the studies presented in this thesis is the determination of the excita- 
tion function of the total cross section for the quasi-free pn — > pnr\' reaction near the 
kinematical threshold. 

Despite the fact that the T| and T|' mesons - which are members of the ground-state 
pseudoscalar nonet [1] - were discovered many decades ago, they are still subject of 
many theoretical and experimantal investigations. 

According to the quark model, T] and r\' mesons can be described as a mixture of the 
singlet and octet states of the SU(3) - flavour pseudoscalar meson nonet. Within the 
one mixing angle scheme, a small mixing angle (0 = — 15.5°) implies that the masses 
of T| and r\' mesons should be almost equal. However, masses of these mesons differ by 
about a factor of two. Additionaly, the mass of the r\' meson does not fit into the SU(3) 
scheme and it is thought to be induced by the gluonic component in its wave function. 
This hypothesis is strenghtened by the decay scheme of mesons like B + or Df since 
the branching ratios of the decay of these mesons into some channels involving r\' are 
significantly higher than the analogous for T| especially in processes requiring the in- 
volvement of gluons [2, 3,4]. 

The properties of the r\' meson should manifest itself in the production mechanism 
in the collisions of nucleons. At present there is not much known about the rela- 
tive contribution of the possible reaction mechanisms of the r\' meson production in 
nucleon-nucleon collisions. It is expected that the r\' meson can be produced through 
heavy meson exchange, through the excitation of an intermediate resonance or via the 
fusion of virtual mesons [5, 6, 7]. However, it is not possible to judge about the mech- 
anism responsible for the r\' meson production only from the total cross section of the 
PP ~^ PP^' reaction [8]. Therefore, one has to investigate the r\' production more de- 
tailed by e.g. selecting separate channels in the relevant degrees of freedom like the 
isospin which means a comparison of both, proton-proton and proton-neutron scatter- 
ing. This conclusion motivated our investigations, which are presented in this thesis. 
A comparison of the close-to-threshold total cross section for the r\' production in both 
PP ~^ PP^i an d pn — > pnr\' reactions constitutes a tool not only to investigate the pro- 
duction of the r\' meson in channels of isospin 1=1 and 1 = but also may provide 
insight into the flavour-singlet (perhaps also into gluonium) content of the r\' meson and 
the relevance of quark-gluon or hadronic degrees of freedom in the creation process. 
It is also possible that the r\' meson is produced from excited glue in the interaction 
region of the colliding nucleons, which couple to the r\' meson directly via its gluonic 
component or through its SU(3)-flavour-singlet admixture [9, 10]. As suggested in ref- 
erence [11], r\' production via the colour- singlet object does not depend on the total 
isospin of the colliding nucleons and should lead to the same production amplitude for 
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the r\ f in the pn — > pnr\' and pp — > ppY reactions. In case of the r\ meson, the ratio 
of the total cross sections for the reactions pn — > pnr\ and pp — > ppr] was determined 
to be T^r, = 6.5 in the excess energy range from ^15 MeV to 160 MeV [12], what 
suggests the dominance of isovector meson exchange in the T| production in nucleon- 
nucleon collisions. Since the fractional amounts of different quark flavours of T| and 
T|' mesons are very similar, in case of the dominant isovector meson exchange - by the 
analogy to the r\ meson production - we can expect that the ratio R^/ should also be 
about 6.5. If however the T|' meson is produced via its flavour-blind gluonium com- 
ponent from the colour-singlet glue excited in the interaction region the ratio should 
approach unity after corrections for the initial and final state interactions. The close- 
to-threshold excitation function for the pp — > ppr\' reaction has been determined in 
previous experiments [13, 14, 15, 16, 17] and the determination of the total cross sec- 
tion for the T|' meson production in the proton-neutron interaction constitutes the main 
motivation for the experiment which is subject of this thesis. 

It is worth to strees that the pn — > pnr\' reaction was never investigated so far. Such 
studies are challenging experimentally because the total cross section is expected to be 
about a factor of fourty to hundred less than in case of the T| meson and additionally 
(in comparison to the r\ near threshold production) the cross section of the multi-pion 
background grows significantly. For some channels (like eg. 3 ft) even by more than 
a few orders of magnitude [18]. The measurement of the pn — > pnr\' reaction is also 
much more difficult in comparison to the pp — > ppr\' reaction due to the neutron in the 
final state and the lack of a pure neutron target. 

The experiment described in this thesis has been performed by the COSY-1 1 col- 
laboration by means of the COSY-1 1 facility at the Cooler Synchrotron COSY at the 
Research Centre Jiilich in Germany. A quasi-free proton-neutron reaction was induced 
by a proton beam impinging on a deuteron target. For the data analysis the proton from 
the deuteron is considered as a spectator which does not interact with the bombarding 
proton, but escapes untouched and hits the detector carrying the Fermi momentum pos- 
sessed at the time of the reaction. The experiment is based on the registration of all 
outgoing nucleons from the pd — > ppnX reaction. Protons moving forward are mea- 
sured in two drift chambers and scintillator detectors and the neutron is registered in 
the neutral particle detector. Protons considered as a spectator are measured by a dedi- 
cated silicon-pad detector. The total energy available for the quasi-free proton-neutron 
reaction can be calculated for each event from the vector of the momenta of the specta- 
tor and beam protons. The absolute momentum of the neutron is determined from the 
time-of-flight between the target and the neutron detector. Application of the missing 
mass technique allows to identify events of the creation of the meson under investiga- 
tion. 

The thesis is divided into ten chapters. The second chapter - following the intro- 
duction - describes briefly the motivation for investigating of the AW — > NNr\' process, 
in particular in view of the study of the T|' production mechanism and its structure. In 
this chapter the most interesting issues concerning the r\' physics are presented, and 
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different possible production mechanisms and the predictions of the total cross section 
for the pn — > pnr\' reaction are given. 

The spectator model is introduced in chapter three, with the description of the quasi- 
free meson production. 

An introduction of the cooler synchrotron COSY and the COSY-1 1 detection setup is 
presented in the fourth chapter. 

The fifth chapter is devoted to the calibration of the detectors. A special emphasis is 
put on the neutral particle and spectator detectors, two new devices which enabled to 
measure the quasi-free pn — > pnX reactions. 

A detailed description of the data analysis and the result of the identification of the 
pn — > pnr\' reaction are presented in chapter six. 

The seventh chapter is devoted to the determination of the luminosity based on the 
quasi-free proton-proton elastic scattering. 

Upper limits of the total cross section as well as of the ratio R^t = c(pn — > pnr\') I 
o(pp — > ppr\') are given in chapter eight. 

The results of the analysis are compared to theoretical predictions in chapter nine. 
The tenth chapter comprises the summary and perspectives. In particular the possibil- 
ity to study of the r|' production in the pure isospin 1 = channel via the pn — > dr\' is 
discussed. 
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2. Motivation 



Understanding of the structure of hadrons is a long standing challenge. Quantum chro- 
modynamics (QCD) is the well established theory of strong interactions that describes 
the underlying forces of coloured quarks which bind them together to form the colour 
neutral hadrons observed in nature. However, due to the increasing coupling constant 
with decreasing energy this theory does not allow at present for the exact description of 
the hadron structure in the low energy domain. One of the hadrons which is especially 
intriguing is the T|' meson. 

The T|' meson was first observed in 1964 [19, 20] in the reaction K p — > At]'. Although 
the T| and r\' mesons - the members of the nonet of the lightest pseudoscalar mesons 
- have been discovered many decades ago, they are still subject of considerable inter- 
est of theoretical as well as experimental studies. Particularly in case of the r\' meson, 
despite more than fourty years of investigations, its structure and properties, as well as 
the production mechanism in collisions involving hadrons are still not well determined. 

According to the quark model, the two physical states of the r\ and r\' mesons are 
considered as a mixture of the SU(3) pseudoscalar octet (rig) and singlet (r]i) states 
with the pseudoscalar mixing angle Q P : 

|r|) =cos9 P |Ti 8 )-sinep|Tii) , 
|iy) =sin0 P |ri 8 )+cos0p|rii) , 

where, following the notation introduced by Gilman and Kaufman [21] and previous 
work by Rosner [22], the SU(3) pseudoscalar octet state T| 8 and singlet state r| i are: 

|t| s ) = —j=\uu + dd — 2ss) , |r]i) = —=\uu + dd + ss) . (2.2) 

The value of the T| — T|' mixing angle in the pseudoscalar meson nonet has been dis- 
cussed many times in the last years [23]. The most up to date value of the mixing angle 
Q P , averaged over all present experimental results, amounts to Qp = — 15.5° ± 1.3° [23]. 
Such small mixing angle implies similar amounts of strange and nonstrange quark con- 
tent: 

|r|) =0.77 --^(uu + dd) -0.63-55 , 

i - (2.3) 

|iy) = 0.63 • ^ (uu + dd) - 0.77 • ss . 

This suggests that the masses of both r\ and r\' mesons should be almost equal. How- 
ever, the values of these masses differ by about a factor of two. Concurrently, the mass 
of the r\' meson does not fit utterly to the SU(3) scheme. More surprisingly, the masses 
of all the pseudescalar mesons, vector mesons and baryons are well described in terms 
of the naive quark model. 
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There are more differences in the physical properties of both mesons. For example: 
excited states of nucleons exists which decay via the emission of the r\ meson, yet none 
of the baryon resonances decay via the emission of the r\' meson [1]. Very high T|' 
apperance in the decays of B and D$ mesons [24, 25, 26] are observed. The branching 
ratio for B + — > K + r\' = (6.5 ± 1 .7) ■ 10~ 5 is much larger than the one for the correspond- 
ing T] channel B + — > K + r\ < 1 .4 ■ 10~ 5 [24]. Similar relations are found in the decay of 

the D 5 meson where ^ p+) = 12.0 ±4.3 exceeds ^f^\ = 4.4±1.2bya 

factor of about three [25, 26]. It is worth to stress that the observed branching ratios do 
not agree with the predictions which ignore the gluonic content of the T|' [2]. 

There exist also essential differences between the production of T] and r\' meson in 
proton-proton collisions close-to-threshold. The total cross section for the pp — > ppr\' 
reaction is by a factor of fourty smaller than the cross section for the pp — > ppr\ reaction 
at corresponding values of the excess energy. The shape of the excitation function is 
also different indicating that the T| meson interaction with nucleons is much stronger 
than the ri'-nucleon one [8]. Thus, it is expected that not only physical properties but 
also the production mechanisms of these mesons should differ from each other. 

As was already discussed, the r|' meson with its mass of m^r ~ 958 MeV is far from 
being "light" and its mass is almost three times larger than the value expected if this 
meson would be a pure Goldstone boson associated with spontaneously broken chiral 
symmetry [27]. The much larger mass of the r\' meson is thought to be induced by the 
non-perturbative gluon dynamics [28, 29, 30] and the axial anomaly Ua(1) [31]. 
A gluonic component of the r\' meson is introduced as a flavour singlet state addition- 
ally to the T]i, which couples directly to the glue [32]. Thus, the x\ f meson can couple 
to gluons not only via the quark and antiquark triangle loop but also directly through 
its gluonic admixture as shown in figure 2.1. 




(a) (b) 

Figure 2.1: Coupling of r\ and r\' to two gluons through quark and anti-quark triangle loop (a) 
and through gluonic admixture (b). Figure and caption are adapted from reference [32]. 

Experimental indications exist that the gluonic content of the r\' meson constitutes 
a significant fraction of the T[' meson wave function. For example, the unexpected 
large branching ratio measured for the decay of beauty particles, B — > r\' +X [24], has 
been interpreted as a possible experimental evidence in this respect [33]. Recenty, the 
KLOE collaboration has estimated the gluonium fractional content of the r\' meson to 
be 14% ± 4% [34] by looking for the radiative decays § — > ri'y and § — > r|y. 
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The most remarkable feature - in the frame of the quark model - distinguishing 
the T|' meson from all other pseudoscalar and vector ground state mesons, is the fact, 
that the T|' is predominantly a flavour-singlet combination of quark-antiquark pairs and 
therefore can mix with purely gluonic states. In any case, if there is a strong coupling 
of the T|' meson to gluons, quark-gluon degrees of freedom may play a significant role 
in the production dynamics of this meson, especially close-to-thereshold where the r\' 
production requires a large momentum transfer between the nucleons and hence can 
occur only at short distances ~ 0.3 fm in the nucleon-nucleon collision. The role of 
gluonic degrees of freedom in the ri'-nucleon system can be investigated for example 
through the flavour- singlet Goldberger-Treiman relation [35]: 

F 2 

2m N G A (0) = Fg Nm >(0) + ^^m^g NNG {0), (2.4) 

which relates the nucleon-nucleon-r]' coupling constant gwr\' with the flavour-singlet 
axial constant Ga(0). The gNNG describes the coupling of the nucleon to the gluons 
arising from contributions violating the Okubo-Zweig-Iizuka rule [36, 37, 38]. The 
coupling constatnt gNNG is in part related to the contribution of gluons to the proton 
spin [39]. As shown by the measurements of the EMC collaboration omission of the 
spin carried by gluons in polarized protons leads to the so-called "spin crisis" [40]. A 
small value of Ga(0) ~ 0.20-0.35 extracted from measurements of the deep inelastic 
muon-proton scattering [40] and a large mass of the r\' meson would be explained by 
the positive gNNG ~ 2.45 value. However, neither gNNr\' nor gNNG have been mea- 
sured directly in experiments so far. There is only an estimation of the upper limit of 
gAWri' < 2.5 derived from the close-to-threshold total cross section for the pp — > ppr\' 
reaction [16]. From this point of view an investigation of processes where the r\' meson 
is produced directly off a nucleon, such as the NN — > NNr\' reaction may be considered 
as a tool for supplying the information about the above coupling constants and the role 
of the gluons in dynamical chiral symmetry breaking in low-energy QCD. 

As suggested by Bass [9, 10] the r\' meson can be created via a contact interaction 
from the glue which is excited in the interaction region of the colliding nucleons. A 
gluon-induced contact interaction contributing to the close-to-threshold NN — > NNr\' 
cross section is derived in the frame of the U(l)-anomaly extended chiral Lagrangian. 
The physical interpretation of the contact interaction is a very short distance (0.2 fm) 
interaction where created gluons could couple to the r\' directly via its gluonic compo- 
nent or through its flavour-singlet admixture T| i . This gluonic contribution to the total 
cross section of the AW — > NNr\' reaction is additional to the production mechanism as- 
sociated with meson exchange [5, 41, 6, 7]. As proposed by Bass in reference [11], the 
T]' meson production via the colour-singlet object does not depend on the total isospin 
of the colliding nucleons (see figure 2.2) and should lead to the same amplitudes of the 
production for the pn — > pnr\' and pp — > ppr\' reactions. This observation motivated the 
studies of the isospin dependence of the r\' meson production in the nucleon-nucleon 
collisions and in particular a measurement of the pn — > pnr\' reaction presented in this 
thesis. 
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Figure 2.2: Diagram depicting pos- 
sible quark-gluon dynamics of the 
AW — > NNr\' reaction. Production via 
a rescattering of a "low energy pseu- 
doscalar pomeron like object" (glu- 
onic current) is isospin independent. 



Treating proton and neutron as different states of the nucleon distinguished only by 
the isospin projection, +^ for the proton and — 5 for the neutron, we may classify the 
AW — > NNX reactions according to the total isospin of the nucleon pair in the initial 
and final state. A total isospin of two nucleons equals 1 for proton-proton and neutron- 
neutron pairs, and may acquire the value 1 or for the neutron-proton system. Since the 
r\' meson is isoscalar, there are only two pertinent transitions for the AW — ► NNX reac- 
tion, provided that it occurs via an isospin conserving interaction. Thus, it is sufficient 
to measure two reaction channels for an unambiguous determination of the isospin 
and 1 cross sections [8]. 



The isospin dependence has been already established in case of the T| meson pro- 
duction, and the total cross sections in both the proton-proton as well as the proton- 
neutron reactions have been measured. In case of the T| meson, the ratio of the to- 
tal cross sections for the reactions pn — ► pnx\ and pp — > ppr\ was determined to be 
i?ri ~ 6.5 [12]. At the lower values of Q, close-to-threshold, the ratio falls down and 
amounts to ~ 3 [42, 43], (see figure 2.3). As explained by Wilkin [44, 45], to a large 
extent this behaviour may plausibly be assigned to the difference in the strength of the 
proton-proton and proton-neutron final state interaction. 
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Figure 2.3: The ratio (7? n ) of the total cross sections for the pn — > pnr\ and pp — » /7/?r( reactions. 
The superimposed line indicates a result of the fit taking into account the final state interaction 
of nucleons [42]. 



Since, 



we have 



G(pp-^ppr\) =0/=i, 
o{pn -> pnTj) = 

O/=0= (2Rn-l)o/=i, 



where / denotes the total isospin in the initial and final state of the nucleon pair. The 
production of the T| meson in this reaction with total isospin 1=0 exceeds the pro- 
duction with isospin / = 1 by a factor of 12. This large difference of the total cross 
sections suggests the dominance of isovector mesons exchange in the creation of the T| 
meson in collisions of nucleons. A recent experimental study of the analysing power of 
the pp — > ppT\ reaction [46] indicates that the n meson exchange between the colliding 
nucleons may be predominant. 

Since the amount of different quark flavours in T| and r\' mesons wave function is 
similar, in case of the dominant isovector meson exchange - by the analogy to the T| 
meson production - we can expect that the ratio ify should be large. If however the 
T|' meson is produced via its flavour-blind gluonium component from the colour- singlet 
glue excited in the interaction region the ratio should approach unity after corrections 
for the initial and final state interactions. The close-to-threshold excitation function 
for the pp — > ppr\' reaction has already been established [13, 14, 15, 16, 17] and the 
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determination of the total cross section for the T]' meson production in the proton- 
neutron interaction motivated the work presented in this dissertation. Recently Bass 
and Thomas [47] argued that the strength of the interaction of r\ and T|' mesons with 
nucleons is sensitive to the singlet-flavour component, and hence to the glue content in 
these mesons. This makes a connection in our endeavour to investigate the structure, 
the production dynamics, and the interaction of the T| and r\' mesons with nucleons [48] . 
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Figure 2.4: Diagrams for the AW — > NNr\' reaction close-to-threshold: 1) - ri'-bremsstrahlung 
(nucleonic current), 2) - "rescattering" term (nucleonic current), 3) - excitation of an intermedi- 
ate resonance (nucleon resonance current), 4) - production through heavy meson exchange, 5) 
- fusion of the virtual mesons (mesonic current). 

Certainly, other production mechanisms which are shown in figure 2.4, such as me- 
son exchange and nucleon resonance currents, must be taken into account before a role 
of gluons in the NNr]' vertex is explored. It is expected that the r\' meson can be pro- 
duced through ri'-bremsstrahlung (nucleonic current), "rescattering" term (nucleonic 
currents), excitation of an intermediate resonance (nucleon resonance current), heavy 
meson exchange or through fusion of the virtual mesons (mesonic curent). The two 
latter mechanisms, which are of short-range, are expected to contribute even more sig- 
nificantly due to the large four-momentum transfer needed between nucleons to create 
the v\' meson. Theoretical studies of the r\' meson production mechanism have shown 
that the existing data could be explained either by mesonic and nucleonic currents or 
by a dominance of two new resonances Sn(1897) and Pn(1986) [5]. Moreover an 
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extended study [41] motivated by the updated data of the yp — > r\'p and pp — > ppr\' 
reactions indicated contributions from resonances Si i( 1650) and Pi i (1870). However, 
it is premature to identify these states, as these authors pointed out. 
Recently, Cao and Lee [7] have studied the near-threshold r\' production mechanism 
in nucleon-nucleon and 71-nucleon collisions under the assumption that the resonance 
Sn(1535) is predominant. In an effective Lagrangian approach which gives a resonable 
description to the pN — > pNr\ and Tip — > pr\ reactions, it was found that the t-channel n 
exchange makes the dominant contribution to the pN — > pNr\' process, and a value of 
6.5 for the ratio of o(pn — > pnr\') to o(pp — > pp"*}') was predicted. 
On the contrary, other authors [6] reproduced the magnitude of the total cross sec- 
tion for the pp — > ppr\' reaction including meson currents and nucleon currents with 
the resonances Sn(1650), Pn(1710) and Pi3(1720). A resonable agreement with the 
pp — > ppr\' data is achieved by the contribution of the meson conversion currents in- 
dicating other terms to be less significant. In the frame of this model the ratio is 
predicted to be 1.5. 

The above considerations shows that our understanding of the r\' meson production 
mechanism is still unsatisfactory. Therefore it is important to test the discussed mecha- 
nisms by confronting them with the experimental results on the isospin dependence of 
the r\' meson production. 
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3. Spectator model 



The absence of pure neutron targets makes measurements of mesons production in 
proton-neutron interactions very difficult, especially very close-to-thereshold where the 
total cross section changes rapidly with the excess energy [49]. Therefore, to measure 
the pn —> pnr\' reaction we used a deuteron target as an alternative. 




Figure 3.1: Schematic picture of the quasi-free pn — > pnr\' reaction. Grey circles represent 
protons, whereas neutrons are depicted as open circles. The Fermi momentum of the nucleons 
inside the deuteron is represented by the dotted arrows. 



Figure 3.1 shows schematically the quasi-free pn — > pnr\' reaction. In the ex- 
periment described in this thesis, a cooled proton beam with a momentum of 
Pbeam = 3.35 GeV/c and a deuteron cluster beam target as a source for neutrons have 
been used. The applied technique is similar to investigations of the n° and T| production 
in proton-neutron collisions carried out by the WASA/PROMICE collaboration at the 
CELSIUS accelerator in Uppsala [12, 50, 51, 52]. 
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Figure 3.2: Total cross section for the pp — > ppr\ reaction as a function of the excess energy 
measured using proton beam colliding with proton target (closed triangles) [53, 54, 55, 56, 57, 
17] and proton beam scattered on deuteron target (open circles) [50]. 

The main difficulty, when using deuterons as a neutron target is the fact that nu- 
cleons inside the deuteron are not at rest but are moving with the Fermi momentum. 
Therefore, in order to achieve the required resolution for the determination of the ex- 
cess energy in the proton-neutron reaction it is mandatory to determine the neutron mo- 
mentum inside the deuteron for each event. The momentum of the neutron is derived 
indirectly based on the measurement of the proton Fermi momentum. For the analysis 
a spectator model is applied. Due to the relatively small binding energy of a deuteron 
(Eg = 2.2 MeV) which is more than three orders of magnitude smaller compared to 
the kinetic energy of the bombarding proton (Tbeam = 2540 MeV), the neutron is con- 
sidered as a free particle in the sense that the matrix element for the quasi-free meson 
production on a bound neutron is equal to the one in free production from an unbound 
neutron. The second assumption of the spectator model used for the evaluation of the 
pn — > pnr\' reaction is that the proton from the deuteron target does not take part in the 
reaction and that it is on its mass shell when the beam hits the target. The validity of 
these assumptions was proven by measurements performed at CELSIUS [50, 51, 52], 
TRIUMF [58] and COSY-TOF [59]. The comparison of the quasi-free and free pro- 
duction cross sections for the pp — > ppr\ reaction done by the WASA/PROMICE col- 
laboration has shown that there is no difference between the total cross section for the 
free and quasi-free process within the statistical errors. The excitation function for 
the pp — > ppf] reaction measured with free and quasi-free proton-proton scattering is 
presented in figure 3.2. Similarly, investigations of pion production carried out at the 
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TRIUMF facility have shown that the experimental spectator momentum distribution 
is consistent with expectations based upon the hypotheses of the spectator model. It 
was also shown, that the magnitude of the differential cross sections for the quasi-free 
pp — > dn + reaction agrees on the few percent level with the free differential cross sec- 
tions. Recently, the validity of the spectator model was proven also by the COSY-TOF 
collaboration. The shape of the angular distribution for the quasi-free np — > pp%~ and 
pn — > pn reactions as well as the form of the momentum distribution of the spectator 
have been measured. The experimental data were consistent with calculations based 
upon the spectator model with an accuracy better than 4% [59]. 




Fermi momentum [MeV/c] kinetic energy [MeV] 

Figure 3.3: Distribution of the momentum (left) and the kinetic energy (right) of nucleons 
inside the deuteron generated according to the analytical parametrisation of the deuteron wave 
function derived from the Paris potential [60, 52]. 

In case of the quasi-free pn — > pnr\' reaction the difference between the on-shell 
and the off-shell total cross section should be even smaller than in case of the quasi-free 
T| meson production since the total energy of colliding nucleons is much larger for near 
threshold T|' creation than for the creation of the r\ meson, whereas the main difference 
between the on-shell and off-shell neutron mass remains the same in both cases. 
The momentum and kinetic energy distribution of a spectator proton is shown in fig- 
ure 3.3. The momentum distribution is peaked at a value of 40 MeV/c but due to the 
long tail the mean momentum is approximately equal to 100 MeV/c. Since the neutron 
momentum changes from event to event, both the total energy in the centre-of-mass 
system (y/s) and the excess energy (Q cm ) vary also and have to be determined for each 
event. 

Figure 3.4 (left) presents the distribution of the excess energy {Qcm) for the ppr\' system 
originating from the pd — > pnr\'p spec reaction calculated with a proton beam momen- 
tum of Pbeam = 3.35 GeV/c and a neutron target momentum smeared out according to 
the Fermi distribution. Here p spec denotes the spectator proton. The broad excess en- 
ergy distribution enables to scan a large range of the excess energy with a single beam 
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momentum setting but in parallel requires the reconstruction of the centre-of-mass en- 
ergy for each event. This can be done only if the four-momentum of target neutrons 
is known, what can be done by measuring the four-momentum vector of the specta- 
tor proton, or by measuring all four-momentum vectors of all outgoing nucleons and 
mesons (proton, neutron and v\'). The first possibility requires the detection of the spec- 
tator proton in a suitable detection unit as it was realized in the COSY-1 1 experimental 
facility. For more details concerning the spectator detector and calibration method see 
chapters 4. and 5. 

As already mentioned, in the frame of the spectator model the proton in the deuteron 
is assumed to be an untouched particle staying on-shell throughout the reaction. Thus, 
by measuring the spectator proton momentum (p spec ) one gets the target neutron mo- 
mentum (p n ) which is equal to —p S pec- The total energy of on-shell proton is equal 
to: 

Espec = \j™ 2 p+P 2 S peci ( 3 • 1 ) 

and hence the total energy of the off- shell neutron can be calculated as: 

E n = m c i E S p ec , (3.2) 



where m p and m c i denote the free proton 
case when the proton is at rest (E spec = 
binding energy of the deuteron. 



and deuteron masses, respectively. Even in 
m p ) the neutron is off mass shell due to the 




excess energy [MeV] mass [MeV/c 2 ] 

Figure 3.4: Left: Distribution of the excess energy (Q cm ) for the pnr\' system originating from 
the pd — > pnr( p S p ec reaction calculated with a proton beam momentum of pteum = 3.35 GeV/c 
and a neutron target momentum smeared out according to the Fermi distribution. Right: Spec- 
trum of the mass of the interacting neutron. An arrow indicates a mass of the free neutron. 

Figure 3.4 (right) shows the distribution of the mass of the interacting neutron as 
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calculated according to the relation: 



m, 



t real 
n 




spec i 



(3.3) 



where E n is the total energy of the off-shell neutron. As it is seen, the max- 
imum of this distribution differs only by 2.2 MeV from the free neutron mass 
(m n = 939.565 MeV/c 2 ). 

The quasi-free meson production is disturbed by some nuclear effects, namely a 
shadow effect and a reabsorption of the produced meson by the spectator proton but 
they are rather of minor importance. In the first case the reduction of the beam flux 
originating in neutron shielding by a spectator proton decreases the total cross section 
by about 4.5% in case of T| meson production [61]. The same effect is expected for the 
r\' meson production. The absorption of the T] meson reduces the total cross section 
by 2% — 4% [61] depending on the energy of the T] meson which was produced. In 
case of a quasi-free T|' meson production this effect is much smaller since the proton-ri' 
interaction is much weaker than the proton-T] one [8]. 
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4. Experimental tools 



The experiment described in this thesis has been performed using the COSY-1 1 facil- 
ity [62, 63, 64] an internal magnetic spectrometer installed at the cooler synchrotron 
and storage ring COSY [65] in Jiilich, Germany. 

4.1 Cooler Synchrotron COSY 

The Cooler Synchrotron COSY stores and accelerates unpolarized as well as polarized 
proton and deuteron beams in the momentum range between 0.6 GeV/c and 3.7 GeV/c. 




Figure 4.1: Scheme of the COoler Synchrotron COSY with indicated internal experiments: 
(COSY-1 1 [62], WASA [66, 67], ANKE [68], PISA [69], EDDA [70]) as well as external 
experiments (JESSICA [71], TOF [72], NESSI [73], GEM [74], BIG KARL [74], MOMO [75], 
HIRES [76]). 



25 



26 



Experimental tools 



Cooling with electrons in the lower momentum range and stochastic cooling in the 
upper momentum regime is used to decrease the spread of the beam momentum and its 
emittance. Proton and deuteron beams are supplied to experiments with the circulating 
beam - "internal experiments" or for experiments with the extracted beam - "external 
experiments". In case of internal experiments, at the highest beam momentum a few 
times 10 10 accelerated particles pass through the target ~ 10 6 times per second. The 
COSY synchrotron - which is schematically depicted in figure 4.1 - has a 180 m cir- 
cumference including two half-circle sections connected by two 40 m straight sections. 
The experiment presented in this thesis is the first ever attempt to measure the cross 
section of the pn — > pnr\' reaction. It has been carried out at the COSY- 11 detection 
setup using a deuteron cluster target and a stochastically cooled proton beam with a mo- 
mentum of 3.35 GeV/c. The experiment was performed in August 2004 and the data 
have been collected for 24 days. It was based on the measurement of four-momenta of 
the outgoing nucleons, whereas an unregistered short lived meson T|' was identified via 
the missing mass technique. 



4.2 COSY-11 facility 



The COSY-11 detector system is schematically depicted in figure 4.2. The detectors 
used for the registration of the spectator proton and the neutron (silicon pad detector and 
neutral particle detector, recpectively) will be described more detailed in the following 
sections. The details of the other detector components used in this experiment and the 
method of measurement can be found in references [62, 64, 77, 78]. Therefore, here 
the used experimental technique will be presented only briefly. 

The COSY-1 1 cluster target, shown in figure 4.3 is installed in front of one of 
the COSY dipole magnets. It can provide streams of clusters of hydrogen (H2) as 
well as of deuteron (D2). The charged products of the measured reactions are bent 
in the magnetic field of the dipole magnet and leave the vacuum chamber through a 
thin exit foil, whereas the beam — due to the much larger momentum — remains 
on its orbit inside the ring. The charged ejectiles are detected in the drift chambers 
(Dl, D2) [62, 64] and the scintillator hodoscopes (SI, S3) [62, 77, 78]. Neutrons and 
gamma quanta are registered in the neutral particle detector [79]. The veto detector 
is used in order to separate neutrons and gamma quanta from charged particles. An 
array of silicon pad detectors (spectator detector) [80] is used for the registration of the 
spectator protons. 
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Figure 4.2: Scheme of the COSY-11 detector system with superimposed tracks from the 
pd — > p sp ec{pnr() reaction. Protons are registered in two drift chambers Dl, D2 and in the 
scintillator hodoscopes SI and S3. The SI detector is built out of 16 vertically arranged mod- 
ules, whereas the S3 hodoscope is a non-segmented scintillator wall viewed by a matrix of 217 
photomultipliers. An array of silicon pad detectors (spectator detector) is used for the registra- 
tion of the spectator protons. Neutrons are registered in the neutral particle detector consisting 
of 24 independent detection units. In order to distinguish neutrons from charged particles a veto 
detector is used. Elastically scattered protons are measured in the scintillator detector S4 and 
position sensitive silicon detector Si mon . Detectors' size and their relative distances are not to 
scale. 
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Figure 4.3: Schematic view of the cluster target. The figure is adapted from reference [77]. 

Protons scattered elastically under large angles are measured in the scintillator 
detector S4 and in another position sensitive silicon pad detector (Si mo „) positioned 
closely behind the S4 counter. 



4.3 Measurement method - general remarks 

The main goal of this experiment is the determination of the total cross section for the 
pn — > pnV[' reaction in the excess energy range between and 24 MeV. In order to 
calculate the total cross section, given in general by the formula: 

C(Qcm) = ~ 777-^ (4.1) 

LxE ef f(Q cm ) 

one needs to determine: i) the efficiency E e /f of the COSY-11 detection system, ii) 
the luminosity L integrated over the time of the experiment, Hi) the excess energy Q cm 
and iv) the number of r\' meson events registered for a given excess energy. In the 
following the experimental conditions, technique of measurement and method of the 
E e ff, L, Qcm, N determination will be presented. 

Events corresponding to the production of at least one charged particle in coinci- 
dence with a neutron (or gamma quantum) were stored on the tape for later analysis. 
Conditions of the main trigger can be written symbolically as: 
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which requires that at least a coincidence of one hit (multiplicity /j > 1) in the neutral 
particle detector (in upper and lower photomultiplier), in the SI detector, and in the 
scintillator wall S3 was demanded. The SI region was restricted to 5 modules (from 
2 nd to 6 th ) in which the proton signal from the pn — > pnr\' reaction was expected. 
During the experiment an additional trigger referred to as T e \ as was set up for the reg- 
istration of the quasi-free p-p elastic scattering. Events which gave a signal in the SI 
detector in coincidence with the S4 detector (see figure 4.2) were accepted as a quasi- 
free proton-proton scattering or elastic proton-deuteron scattering. 
For each charged particle, which gave signals in the drift chambers, the momentum vec- 
tor at the reaction point can be determined. For the analysis of the pn — > pnr\' reaction, 
first the trajectories of the outgoing protons are reconstructed [81] based on signals in 
the drift chambers Dl and D2 - and then knowing the magnetic field of the dipole mag- 
net - its momentum vector is determined. Independently, the velocity of the particle 
is derived based on the time-of-flight measured between the SI and the S3 detectors. 
Knowing the velocity and the momentum of the particle, its mass can be calculated, 
and hence the particle can be identified. After the particle identification the time of 
the reaction at the target is obtained from the known trajectory, velocity, and the time 
measured by the SI detector. 

The neutral particle detector delivers the information about the time at which the 
registered neutron or gamma quantum induced a hadronic or electromagnetic reaction 
in the detector volume, respectively. The time of the reaction combined with this in- 
formation allows to calculate the time-of-flight (TOF^) of the neutron (or gamma) 
between the target and the neutral particle detector, and — in case of neutrons — to 
determine the absolute value of the momentum (p), expressed as: 

/ 1 

where m n denotes the mass of the neutron, and / stands for the distance between the 
target and the neutral particle detector. The direction of the neutron momentum vector 
is deduced from the angle defined by the centre of the hit module inside the neutral 
particle detector. 

Similarly, to determine the direction of the momentum vector of the spectator pro- 
ton, the angle defined by the centre of the hit segment inside the spectator detector is 
used. In order to calculate the four-momentum of the spectator proton its kinetic en- 
ergy (T spec ) is directly measured as the energy loss in the silicon detector. Knowing the 
proton kinetic energy (T spec ) one can calculate its momentum using the relationship: 



p=\J(T spec + m p ) 2 -m p 2 , (4.3) 

where m p denotes the proton mass. 

To evaluate the luminosity, the quasi-elastic proton-proton scattering is measured 
during the data taking of the pn — > pnr\' reaction, with one proton detected in the drift 
chambers and the scintillator hodoscopes and the other proton registered in the silicon 
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detector Si mon . The elastically scattered protons can be well separated from the multi- 
particle reactions, due to the two body kinematics. 

Based on Monte Carlo studies, the beam momentum value, the position of the spec- 
tator detector and the configuration of the neutral particle detector have been optimized 
before the experiment in order to maximize the detection efficiency, the resolution of 
the excess energy determination and in order to achieve a relatively high missing mass 
resolution [82]. 

In the next two chapters the method of the detector calibration and evaluation of 
the data, namely determination of the excess energy Q cm , number of T|' mesons cre- 
ated in proton-neutron collisions, the luminosity L and the detection efficiency will be 
presented. 



5. Calibration of the detectors 



5.1 Time-space calibration of the drift chambers 

In this section we give an account on the time-space calibration of the drift chambers. 
The drift chamber Dl (see figure 4.2) contains six detection planes, two planes have 
vertical wires and four have inclined wires. The second drift chamber (D2) contains 
two more planes with vertical wires (altogether eight). The wires in consecutive pairs 
of planes are shifted by half of the cell width in order to resolve the left-right position 
ambiguity with respect to the sense wire. Both chambers were operating with 50%- 
50% argon-ethane gas mixture at atmospheric pressure. 



0.05 




drift time + offset [ns] drift time + offset [ns] 

Figure 5.1: Left: Example of the spectrum used for the drift chamber calibration. The shaded 
area presents mean values of AX as a function of the drift time. Upper and lower histograms in- 
dicate a band of one standard deviation of the AX distribution. The meaning of AX is explained 
in the text. Right: Distance from the sense wire as a function of the drift time for an arbitrarily 
chosen plane of the Dl chamber. 

When a charged particle passes through the drift cell it invokes an ionization of 
the gas mixture. The positively charged ions drift to the cathode wires with negative 
potential, whereas the free electrons drift towards the sense wire with positive potential. 
The drift time of electrons is measured and hence the distance between the particles 
trajectory and the sense wire can be determined. A relation between the drift time and 
the distance from a sense wire is established from experimental data. Due to variations 
of the drift velocity caused by changes of the pressure, temperature and humidity of the 
air inside the COSY tunnel, the drift chamber calibration has to be performed for time 
intervals not longer than few hours. The calibration is derived in an iterative way. First, 
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having an approximate time-space function x(t) for a given sample of measured events 1 , 
a distance of particle trajectory from a sense wire is calculated, and a straight line is 
fitted to the obtained points. Next, the deviation AX between the fitted and measured 
distances of the particle's trajectory to the wire is calculated. An example of a mean 
value of AX as a function of the drift time is shown in figure 5.1 (left). The shaded 
area represents the mean values of AX. Upper and lower histograms indicate a band 
of one standard deviation of the AX distribution. Afterwards the relation between the 
drift time and distance from the wire is corrected by the determined mean value of AX. 
The procedure is repeated until corrections become negligible compared with expected 
position resolution of the chambers. 

Figure 5.1 (right) shows the time-space calibration of an arbitrarily chosen plane of the 
Dl chamber. The nearly linear dependence is seen for the drift time range from 100 ns 
up to 500 ns. 

5.2 Time calibration of the SI and S3 hodoscopes 

In this section a method of the time-of-flight calibration is presented. 

The time-of-flight measurement on the known distance between the SI and the 
S3 hodoscopes enables to calculate the velocity of the particle which crosses both de- 
tectors. In order to calculate the time-of-flight, relative time offsets between each 
detection unit of the SI and S3 detectors have to be established. 

The SI detector (see figure 4.2) is built out of sixteen scintillation modules with 
dimensions 45 cm x 10 cm x 0.4 cm. The modules are arranged vertically with a 
1 mm overlap and are read out at both ends by pairs of photomultipliers. 
The start time for the time-of-flight measurement is calculated as the mean value of 
times measured by the upper and lower photomultiplier of the S 1 detector module. The 
measured TDC values for a single SI module may be expressed as: 

TDC U S [ = t S i + ^-+ tZf + off set up ~ T trigger , 

PL 

TDC d s r = t S i + ^+tZ lk + offset dw - T trigger . (5.1) 

PL 

The tsi is the real time when a particle crosses the detector, y denotes the distance 
between the upper edge of the active part of the detector and the hit position, t^ lk and 

tdw k stan d for the time walk effect, and offset up , offsetdw comprise all delays due to 
the utilized electronic circuits. The common start signal for all TDC modules is denoted 
by the T tr i gger and it is the same for SI and for S3 detector. Thus the mean value of the 

! In case of this experiment the first approximation was taken from the previous COSY-1 1 runs. It 
is worth noting that in general an approximate drift time to drift distance relation can be determined by 
integration of the drift time spectra as provided by the uniform irradiation method [83]. 
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time measured by the SI module is given by: 

TDC U J\ + TDCi? 

tsi = SJ -^ — " offsets, - t w s f k + T triggen (5.2) 

where offsets, comprises all constatnt terms from eq. 5.1. 

The S3 hodoscope is built out of a non-segmented scintillator wall with the dimensions 
220 cm x 100 cm x 5 cm viewed by a matrix of 217 photomultipliers. The stop signal 
for the time-of-flight measurement is calculated as the average of times obtained from 
all the photomultipliers that produced a signal. The TDC value for a single photomul- 
tiplier is expressed similarly as in the case of the SI detector: 

TDC S 3 = t S 3 + C + t^ lk + offsets^ ~ T trigger , (5.3) 

where ts3 denotes the real time when a particle crosses the detector and the t^ s is 
the time needed by the light signal to pass from the scintillation origin down to the 
photomultiplier photocathodes. Applaying equations 5.1 and 5.3 one can calculate the 
time-of-flight between the SI and S3 hodoscopes as: 

TOF s3 ~ sl =ts3-tsu 

T T)C U P _i_ T ptfdw 

Tn rS3-Sl Tnr pos walk n ff„„ t SI 51 , +walk , „ ff „„+ 

LOb = 1 DCs3 - t S3 - t s3 - offsets? ^ ^ f Sl + °JJ set Sl ■ 

(5.4) 

Since in both equations (5.1, 5.3) the T tr i gger time is the same, the TOF s ~ 1 51 is inde- 
pendent of the triggering time. The tj£ s as well as t^f k and t^ lk may be calculated 
based on the known calibration parameters and amplitude of the signals [77]. Thus the 
only unknown quantities in equation 5.4 for TOF s ~ 1 are the time offsets for the indi- 
vidual photomultipliers of S3 detector (offsets?,) and time offsets for a single module 
of SI detector (offsetsi). 

In order to establish these effects we have selected events with only one track re- 
constructed in the drift chambers, with signals in only one SI module and signals in 
the S3 hodoscope for those photomultipiers which are at the position expected from the 
extrapolation of the particle trajectory determined in the drift chambers. 
In the first approximation the relative offsets of the S 1 modules were established com- 
paring the time between neighbouring modules for signals from particles which crossed 
the detector through the overlapping regions. Assuming, that the time offsets for the SI 
detector are approximately correct, we have established the time offsets for each pho- 
tomultiplier of S3 detector, by comparing the time-of-flight values calculated from the 
time signals in the SI and S3 detectors (see eq. 5.4) with the time-of-flight determined 
from the reconstructed momenta of the particles in the magnetic field. Having the time 
offsets for the S3 hodoscope adjusted, we have again re-calculated the offsets for the 
SI detector in the same way as it was done for the S3 detector. This procedure was 
repeated a few times until the corrections of the offsets became negligible. 
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Figure 5.2: The difference between the time-of-flight (TOFsisi) determined from signals 
registered in the SI and S3 scintillators and the time-of-flight calculated from the reconstructed 
momenta of particles versus the photomultiplier ID of S3. The figure shows result after the S3 
hodoscope calibration. Photomultipliers which are positioned at the lower and upper edge of 
the scintillation wall have registered essentially less events (ID = 130 and 139). 

Figure 5.2 presents the distribution of difference between the time-of-flight 
(TOFs3si) on the S1-S3 distance determined from signals registered in SI and S3 
scintillators and the time-of-flight calculated from the reconstructed momenta of par- 
ticles for photomultipliers no. 130-149 as obtained after the third iteration. 



5.3 The neutral particle detector 

The installation of the neutral particle detector at the COSY-1 1 facility enabled to study 
reactions with neutrons in the exit channel. It allows, for example to investigate quasi- 
free meson production in proton-neutron interactions or charged hyperon production 
like E + via the pp — > nK + Z + reaction [84, 85]. This detector [86] delivers the time 
at which neutron or gamma quantum induces a hadronic or electromagnetic reaction 
inside the detector volume, respectively. This information combined with the time of 
the reaction at target place — deduced using other detectors — enables to calculate the 
time-of-flight between the target and the neutral particle detector and to determine the 
absolute momentum of the registered neutrons. 

In this section a method of time calibration will be demonstrated and results achieved 
by its application will be presented and discussed. 
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Information about the deposited energy is not used in the data analysis because the 
smearing of the neutron energy determined by this manner is by more than an order of 
magnitude larger than established from the time-of-flight method. The experimental 
precision of the missing mass determination of the pn — > pnr\' reaction [87, 88] strongly 
relies on the accuracy of the reconstruction of the momentum of neutrons, therefore the 
time calibration of the neutral particle detector has to be done with high precision. 

Previously, the neutron detector consisted of 12 detection units [89], with light 
guides and photomultipliers mounted on one side of the module. In order to improve 
the time resolution of the detector additional light guides and photomulipliers were in- 
stalled, such that the light signals from the scintillation layers are read out at both sides 
of the module. The neutral particle detector used for the experiment described in this 
thesis consists of 24 modules, as shown in figure 5.3. Each module is built out of eleven 
plates of scintillator material with dimensions 240 mm x 90 mm x 4 mm interlaced 
with eleven plates of lead with the same dimensions. The scintillators are connected at 
both ends of a module to light guides — made of plexiglass — whose shape changes 
from rectangular to cylindrical, in order to accumulate the produced light on the circular 
photocathode of a photomultiplier. There the light pulse is converted into an electrical 
signal, which is provided to the ADC and to the TDC converters. 
The neutron detector is positioned at a distance of 7.36 m from the target in the config- 
uration schematically depicted in figure 5.3. The detector covers the laboratory angular 
range of ±1.84° in x and ±1.1° in y direction. 



Figure 5.3: Left: Scheme of the single detection unit of the COSY-11 neutral particle detector. 
LG and PM denotes light guide and photomultiplier, respectively. The figure shows also the 
definition of variables y and L used in the text. Right: Configuration of the detection units of 
the neutral particle detector. 

This configuration and the distance to the target has been chosen to optimize the 
efficiency for the reconstruction of the pn — > pnr\' reaction [82]. The choice was based 
on results from Monte Carlo studies of the acceptance and efficiency as a function of the 
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distance and configuration [82]. The choice for the thickness of the scintillator plates 
was also based on the simulation studies [90] . 
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Figure 5.4: Efficiency of the 20 cm thick neutron detector as a function of neutron kinetic 
energy simulated for various layer thicknesses of the scintillator and iron plates. The solid line 
with triangles is for pure scintillator material. Diamonds denote a layer thickness of 0.5 mm, 
closed circles of 2 mm, open circles of 5 mm, and crosses of 25 mm. The figure has been 
adapted from ref. [90]. Note that for better visualisation the vertical axis starts at 20%. 

As can be deduced from figure 5.4 the maximum efficiency, for a given total thick- 
ness, for the registration of neutrons — in the kinetic energy range of interest for the 
pd — > pnr\'p sp reaction (~ 300 MeV - ~ 700 MeV) — would be achieved for a homo- 
geneous mixture of lead and scintillator. However, in order to optimize the efficiency 
and the cost of the detector the plate thickness has been chosen to be 4 mm. This results 
in an efficiency which is only of few per cent smaller than the maximum possible. The 
functioning of the detector was already tested in previous experiments [42, 79, 85]. 

5.3.1 Time signals from a single detection unit 

A schematic view of a single detection module of the COSY-1 1 neutral particle detector 
is shown in figure 5.3 (left). The time (T exp ) from a single module is calculated as an 
average time measured by the upper and lower photomultiplier. Namely: 



TDC up + TDC 
~ 2 



(5.5) 
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where TDC up ' denotes the difference between the time of signals arrival from the 
photomultiplier and from the trigger unit to the Time-to-Digital-Converter (TDC). 
This can be expressed as 2 : 

TDC up = t real + offset"* + — - Tagger, (5.6) 



TDC dw = t real + offset^ + - T tri gg er , (5.7) 

where L stands for the length of a single module, y denotes the distance between the 
upper edge of the active part of the detector and the point at which a neutron induced the 
hadronic reaction, t real is the time at which the scintillator light was produced, T tr i gger 
represents the time at which the trigger signal arrives at the TDC converter, and cl 
denotes the velocity of the light signal propagation inside the scintillator plates. The 
parameters offset" p and offset dw denote the time of propagation of signals from the 
upper and lower edge of the scintillator to the TDC unit. 

Applying equations 5.5, 5.6 and 5.7 one can calculate a relation between T exp and 

^real. 

. offset"* + offset^ + 

T e Xp = f eal + _ £L _ = f eal + offset _ (5 g) 

where the value of "offset" comprises all delays due to the utilized electronic circuits, 
and it needs to be established separately for each segment. It is worth noting, that due 
to the readout at both ends of the detector the T exp is independent of the hit position 
along the module, as it can be deduced from equation 5.8. 



5.3.2 Relative timing between modules 

Instead of determining the values of "offsets" from equation 5.8 for each detection unit 
separately, the relative timing between modules will be first established and then the 
general time offset connecting the timing of all segments with the S 1 detector will be 
found. 

In order to establish relative time offsets for all single detection units, distributions of 
time differences between neighbouring modules were derived from experimental data 
and compared with simulated distributions. A time difference measured between two 
modules can be expressed as: 

Ay = T- xp - T- xp = t™ al - t] eal + (offset; - offset;) , (5.9) 

where T[ xp and T exp stand for the time registered by the i th and j th module, respectively. 
Examples of A ;J - spectra determined before the calibration are presented in figure 5.5. 



2 Here we omit the time walk effect which will be discussed later in section 5.3.3 
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Figure 5.5: Distributions of the time difference between the 8 and the 13 , the 9 and the 
14 th , and the 18'' 1 and the 23 th module of the neutron detector, as determined before the calibra- 
tion assuming that offsets are equal to zero. 



One can note that the peaks are shifted from the zero value and additionally the dis- 
tributions contain tails. The tails reflect the velocity distribution of the secondary par- 
ticles. In order to determine reference spectra corresponding time differences between 
the modules were simulated using the GEANT-3 code. The result of the simulation is 
shown in figure 5.6. To produce these spectra the quasi-free pn — > pnr\' reaction has 
been simulated. 




Figure 5.6: Simulated distributions of the time difference between the 8 and the 13 , the 9' 
and the 14 th , and the 18 th and the \3 th module of the neutron detector. 



The values of the relative time offsets were determined using a dedicated program 
written in Fortran 90 [91, 92]. It adjusts values of offsets such that the time difference 
obtained from the experimental data and from simulations equals to each other for each 
pair of detection units. Furthermore, from the width of the spectra one can obtain the 
information about the time resolution of a single module, which was extracted to be 
0.4 ns [91]. 
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Figure 5.7: Distributions of the time difference between the 8 and the 13 , the 9 and the 
14 th , and the \8 th and the 23 th module of the neutron detector, as obtained after the calibration. 

Figure 5.7 presents experimental distributions of time differences between neigh- 
bouring modules as determined after the calibration. Examining figures 5.7 and 5.6 it 
is evident that now the peaks are positioned as expected from simulation. 

5.3.3 General time offset and the time walk correction 

The momenta of neutrons are determined on the basis of the time between the reaction 
and the hit time in the neutral particle detector. The time of the reaction can be deduced 
from the time when the exit proton from the pn — > pnr\' reaction crosses the S 1 detector 
(see fig. 4.2) because the trajectory and the velocity of this proton can be reconstructed. 
Therefore, to perform the calculation of the time-of-fiight between the target and the 
neutral particle detector, a general time offset of this detector with respect to the S 1 
counter has to be established. This can be done by measuring the gamma quanta for 
which the velocity is constant. 

The time-of-fiight between the target and the neutral particle detector (TOF N ) 
is calculated as the difference between the time of the module in the neutral particle 
detector which fired as the first one (t^ al ) and the time of the reaction in the target 
(t r r eal ): 

TOF N = t r N eal -t r r eal (5.10) 

The time of the reaction is obtained from backtracking the protons through the known 
magnetic field and the time measured by the SI detector (tsi). Thus t r r eal can be ex- 
pressed as: 

f r eal = tsi - TOF sl = TH P - offset 81 + T trigger - TOF sl , (5.11) 

where TOF Si denotes the time-of-fiight between the target and the SI counter, 
offset 81 denotes all delays of signals from the SI detector, and analogically as for 
equation 5.8 the ts\ and T^ p denote the time at which the scintillator light in the SI 
detector was produced and the averaged time registered by the TDC unit (see eq. 5.5), 
respectively. Since the time in the neutron detector reads: 

f eal = T^ p -off S e^ + T trigger (5.12) 
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we have: 

TOF N = T™ p + TOF sl -T™ p -offset N +offset sl = T N xp + TOF sl -T** p + offset G 

(5.13) 

By offset the general time offset of the neutron detector with respect to S 1 is denoted. 
In order to determine the value of offset the pd — > A + B + yX reaction is used, where y 
is measured by the neutral particle detector and for which the TOF N between the target 
and the neutral particle detector is known. A + and B + represent any two positively 
charged particles (7t + 7t + , pn + , pd, dn + ). Events corresponding to the pd ^A + B + yX 
reaction have been identified by measuring the outgoing charged as well as neutral 
ejectiles. Positively charged particles (pions, protons and deuterons) were detected by 
means of the drift chambers (Dl, D2) and scintillator hodoscopes (SI, S3). Gamma 
quanta originated predominately from the 7t° meson decay are registered in the neu- 
tral particle detector. Knowing the time-of-flight for gamma quanta Ty on the distance 
between the target and the neutral particle detector the general time offset was deter- 
mined by comparing the measured time with the nominal value of Ty. The measured 
time-of-flight spectrum is shown in figure 5.8. Ty should be equal to 24.5 ns. As can 
be inspected from the figure 5.8, the time depends on the amplitude. This phenomenon 
is occuring for measurements using a leading-edge discriminator and is referred to as 
the time- walk effect [93]. In the case of a constant threshold value, two signals of 
different pulse hights but exactly coincident in time may trigger the discriminator at 
different times. An offline correction can be applied to minimize this effect assuming 
a linear dependence between the time walk and the inverse of the square root of the 
signal charge [94]. The function used for offline time- walk correction was 

t c = t m -a + $( 7=j= + / 1 i )» ( 5 - 14 ) 
VADC u p VADC dw 

where t c and t m are corrected and measured times, and a and (3 are coefficients deter- 
mined from the data, and ADC up and ADC dw denote the charge of the signal measured 
at the upper and the lower end of the detection module. 

Figure 5.8 shows the dependence of the time-of-flight for neutral particles as a 
function of values of ADC signals in the neutral particle detector. The time-of-flight 
shown in the figure was measured between the target and the neutral particle detector 
before (left) and after (right) the time walk correction. The spectrum was obtained 
under the condition that in coincidence with a signal in the neutral particle two charged 
particles were registered in the drift chambers. The time in the neutral particle detector 
was corrected in such a way as if the interaction point of particles was in the first row 
of the detector modules only. A clear signal originating from the gamma rays is seen 
over a broad enhancement from neutrons. This histogram shows that a discrimination 
between the signals originating from neutrons and gamma quanta can be done by a cut 
on the time of flight. For the gamma quanta the time-of-flight value is independent of 
their energy. This fact was used to determine the a and (3 coefficients. The right panel 
of figure 5.8 presents the analogous distribution for which a corrected time t c was used 
for the time-of-flight calculation instead of the t m . 
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Figure 5.9: Distribution of the time-of-flight between the target and the neutral particle detec- 
tor before (left) and after (right) time walk correction. 

Figure 5.9 presents the projection of the distributions from figure 5.8 on the time- 
of-flight axis. Both distributions of the time-of-flight were fitted by a Gaussian func- 
tion in the range where the clear peak originating from gamma quanta appears. The area 
around the peak was fitted by a second order polynomial. The left panel of figure 5.9 
presents the distribution of the time-of-flight before the time walk correction, and the 
fit was resulting in a time-of-flight resolution of a = 1.2 ns. The time-of-flight for 
which the corrected time in the neutral particle detector was used is shown in the right 
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panel of figure 5.9. The application of the time walk correction improved the time res- 
olution to a = 0.6 ns. It is worth noting that this is an overall time-of-flight resolution 
resulting from the time resolution of the neutron and SI detectors and the accuracy of 
the momentum reconstruction of charged particles, needed for the determination of the 
time of the reaction in the target. 

Figure 5.10 (left) shows the TOF distribution as a function of modules number. As 
expected from the known absorption coefficients [95], the gamma quanta are predom- 
inantly registered in the first row of the detector (see fig. 5.3) whereas the interaction 
points of the neutrons are distributed more homogeneously. Indeed, a clear signal orig- 
inating from the gamma quanta is seen on the time-of-flight distribution, when taking 
into account signals from the first row of the neutron detector only, as it is shown in 
figure 5.10 (right). 




10 20 30 40 50 60 70 80 

time of flight [ns] 



350 : 
300 : 
250 : 



i2 200 



100 
50 




10 20 30 40 50 60 

time of flight [ns] 



Figure 5.10: Left: Time-of-flight determined between the target and the neutron detector as a 
function of modules's number. Right: TOF N distribution determined taking signals from first 
row of the neutral particle detector only. 



Therefore, in order to raise the confidence to the determination of the general offset 
we have taken signals in the first row of the neutral particle detector. The data have 
been restricted only for the pd — > pdyX channel. In this case due to the baryon number 
conservation there is only one possible source of a signal in the neutral particle detec- 
tor, namely a gamma quantum, which originates predominantly from the decay of n° 
mesons in the target. 
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Figure 5.11: Time-of-flight distribution determined between the target and the neutron detector 
obtained under the assumption that additionally to a signal in the neutron detector one proton 
and one deuteron were identified from the signals in the drift chambers. 

Figure 5.11 presents the experimental distribution of the time-of-flight between the 
target and the neutron detector for neutral particles, with the requirement that addition- 
ally to a signal in the neutron detector two charged particles were registered and that 
one of them was identified as a proton and the other as a deuteron 3 . As expected, in this 
spectrum only a signal from the gamma quanta is seen. The clear peak is positioned at 
the value equal to 24.5 ns, which corresponds to the time-of-flight for the light on the 
target-neutral particle detector distance. 



5.3.4 Efficiency determination 

The efficiency of the COSY-1 1 neutral particle detector - which is an important factor 
for determining the absolute values of cross sections - was determined using two in- 
dependend simulation programs. For the first case, a procedure based on the GEANT- 
3 (GEometry ANd Tracking) code [96] was used for the simulation of the hadronic 
cascades induced in matter by neutrons. The same procedure was repeated using the 
FLUKA 4 (FLUktuierende KAskade) [97, 98] simulation program. 

The efficiency of the neutron detector is given by the ratio of the number of events, 
for which an energy deposited in the scintillator material was larger than the threshold 

3 The method of identifying the particles registered by the drift chambers and scintillator hodoscopes 
is described in Chapter 6. 

4 The simulations were performed with the 2008 version. 
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value at least in one of 24 detection units to the number of generated neutrons. The 
value of the calculated efficiency as a function of the kinetic energy of the neutrons 
is shown in fig. 5.12 (left). Open squares denote result obtained using the GEANT-3 
package and the outcome of the simulation using FLUKA-2008 is presented as black 
circles [99]. 
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Figure 5.12: Left: The efficiency distribution as a function of the kinetic energy of the neutrons 
determined assuming that the threshold is equal to MeV. Right: The relation between the 
threshold value and the efficiency for neutrons with an energy of 300 MeV. 

The kinetic energy of the neutrons from the pn — > pnx( reaction varies from 
300 MeV up to 700 MeV for the 3.35 GeV/c beam momentum, and as can be in- 
spected from fig. 5.12 (left) the efficiency is fairly constant in this range. It is important 
to mention that the results determined with the different programs are in agreement 
within ±3%, independently of the threshold value. 

Also studies of the variation of the efficiency depending on the threshold have been 
conducted. In the experiment the threshold was set to about 0.1 MeV and therefore the 
values from up to 0.6 MeV were scanned. The result is presented in fig. 5.12 (right). 
For both, the GEANT and FLUKA-2008 simulation the efficiency changes by about 
10% over the 0.6 MeV range of the threshold. 



5.4 Veto detector 

In order to distinguish between charged and neutral particles, in front of the first layer 
of the neutral particle detector an additional scintillator detector was installed as it is 
shown in figure 5.13. This detector - referred to as veto detector - is built out of four 
overlapping modules with dimensions of 400 mm x 200 mm x 4 mm. The lightguides 
and photomultipiers are mounted at the upper and lower edge of the modules such that 
the light signals are read out at both sides of the module. 
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It permits to reject the background originating from charged particles hitting the neutral 
particle detector. The discrimination between signals arising from the only possible re- 
maining neutrons and gamma quanta is done by a cut on the time-of-flight (see chapter 
6.). 

VETO 




x 



z 



Figure 5.13: Layout of veto detector and neutral particle detector. 




Figure 5.14: Left: The relation between the time-of-flight between the target and the neutral 
particle detector (TOF N ) and the time-of-flight between the target and veto detector (TOF v ). 
Right: ADC signals in the scintillator veto detector. The red arrow shows the ADC range for 
charged particles, which are rejected in the further analysis, whereas the blue arrow point to the 
ADC range for which events are interpreted as neutral particles. Noise signals are extinguished 
by the coincidence between the upper and lower photomultiplier. 

As a first step of the time calibration of the veto detector the time-of-flight between 
the target and the veto detector has been calculated for each module separately as sum- 
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ing that all time offsets are equal to zero. Next, the values of the calibration constants 
(offsets) were adjusted such that the time-of-fiight value is equal for each module. 
Figure 5.14 (left) shows the relation of the time-of-flight between the target and the 
neutral particle detector (TOF N ) and the time-of-flight between the target and the veto 
detector (TOF v ). In principle this relation should be linear, however we observe events 
beyond the correlation line with time signals in the veto detector larger than in the neu- 
tral particle detector. This may happened when a neutral particle (which does not give a 
signal in the veto detector) induces electromagnetic or hadronic reactions in the neutral 
particle detector, and charged products of these reactions scattered backwards resulting 
in signals in the veto detector. This effect is reflected in figure 5.14 (right) which shows 
the ADC signals in one of the veto detector module. Two well separated peaks are 
observed. The peak in the range from 30 to 300 channels originates from events with 
charged particles which crossed the veto detector and then entered the neutral particle 
detector. When taking into account events which are above the correlation line only, 
the peak in the range between channel numbers 30 and 300 vanishes and only events 
with small ADC value up to 30 remain. Therefore, to discriminate between charged 
and neutral particles a cut on the veto detector ADC was applied, which is denoted in 
the figure 5.14 (right) by a dotted line. 



5.5 Spectator detector 

In order to identify the unobserved r\' meson in the quasi-free pn — > pnX process by 
means of the missing mass technique, one has to measure not only the momenta of out- 
going proton and neutron, but also a determination of the pn reaction energy (y/s) for all 
measured events is needed. The latter was realised by measuring the four-momentum 
vector of the spectator proton. 



5.5.1 Scheme of the spectator detector 

Before installation at the COSY- 1 1 detection setup, the spectator detector has been used 
by the PROMICE/WASA collaboration at the CELCIUS storage ring of The Svedberg 
Laboratory [80]. Detectors which are located inside the scattering chamber, close to 
the interaction point, must be compatible with high vacuum and resistant to radiation 
environment. Moreover, a very good energy resolution is highly required in order to 
determine the four-momentum of the spectator protons. These factors exclude the use 
of plastic scintillator. Solid state devices are well suited for this purpose, therefore sili- 
con was chosen as an active material in the spectator detector. 
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Figure 5.15: Schematic view of the spectator detector's single modules. Left: front layer (dE). 
Right: back layer (E). Indicated dimensions are given in mm. 

The spectator detector consists of four double-layered modules. The layers are 
produced of 300 thick silicon. Each of the front layers, which are the planes closer 
to the beam, contain eighteen silicon pads. The active area of a single pad amounts 
to 20 mm x 5 mm. The back layer contains six silicon pads, with an active area of 
20 mm x 18 mm. The schematic picture of the front (dE) and the back (E) layers is 
shown in figure 5.15. 




Figure 5.16: Schematic view of the spectator detector as installed at the COSY- 11 facility. The 
z-axis corresponds to the direction of the beam, x is chosen parallel to the acceleration plane 
and y is vertical. 

The silicon detector system is placed in the high vacuum at the COSY storage ring 
at a distance of about 5 cm upstream from the interaction point as it is schematically 
shown in figure 5.16. The active area covers about 22% of the full solid angle. The 
arrangement was a compromise between the technical needs for the installation and the 
angular resolution. 
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5.5.2 Energy calibration and position optimization 

During the experiment, the charge of the signal induced in the detector by the particles, 
which is proportional to the energy deposited in the detection unit, has been measured. 
Denoting by ADC \ and ADC2 the measured charge in the first and second layer, respec- 
tively, the real energy loss can be expressed as: A E\ = a\ADCi and A E2 = QL2ADC2, 
where oci and 0C2 are the calibration constants for two layers. The aim of the spectator 
detector calibration was to find the exact position of the detector and then to establish 
the calibration constants for each of the 96 silicon pads. The position of the spectator 
detector inside the beam pipe is described by only three parameters [82]. Variable dz 
denotes the translation along the z-axis in the main reference system (MARS). Param- 
eters dl and alpha account for translation and rotation of the flange used for supporting 
the spectator detector. 

Differences in shapes of the dE\(dE2) distributions are due to the fact that the ef- 
fective thickness of the silicon modules of the detector "seen" by particles outgoing 
from the interaction point depends on the particles' incident angle on the detector sur- 
face. These variations allowed to find the exact position and orientation of the specta- 
tor detector inside the COSY-ring. Figure 5.17 presents an example of the calculated 
dE\ {dE%) distributions for a few arbitrarily chosen sets of parameters dz, dl, and alpha. 
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Figure 5.17: Energy losses in the first layer versus energy losses in the second layer as calcu- 
lated for different positions of the spectator detector. The values of dz, dl are given in cm and 
alpha is given in rad. 



Assuming that the set of calibration constants (oci, OC2) are correct, experimental 
distributions of the dE\{dE2) have been fitted to the calculated energy losses, and ex- 
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amined by varying the parameters describing the position of the spectator detector. Af- 
ter the determination of the new position of the detector the experimental data points on 
the dE\ (dE2) plot were fitted to the expected dE\ {dEj) function, with calibration con- 
stans oci and 0C2 as free parameters for each detection pair. The procedure was repeated 
until the changes become negligible. 




Figure 5.18: Energy losses in the first layer versus the second layer as measured at COSY-11 
with a deuteron target and a proton beam with a momentum of 3.35 GeV/c before (left) and 
after (right) the calibration. Points represent the experimental data, whereas the theoretical 
calculations are indicated by the solid curve. 



Figure 5.18 (left) shows the energy losses in the first layer versus the second layer 
before the calibration constants oci and 0C2 were determined. The theoretical calcula- 
tions are indicated by the solid curve. Fig. 5.18 (right) shows the corresponding plot 
after the calibration. The position of the spectator detector has been determined with 
an accuracy of ±1 mm, consistent with the nominal values based on the geometrical 
designs of the setup and with the values obtained for the analysis of the pn — ► pnr\ test 
reaction [42]. Taking into account the size of the stream of the deuteron target with a 
diameter of 9 mm [100] the accuracy of the position determination within ±1 mm is 
satisfactory. 

Figure 5.19 shows the distribution of the excess energy with respect to the pnr\' system 
determined for the pn reaction. The solid histogram was obtained with the position of 
the spectator detector determined in the above described procedure. 
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Figure 5.19: The distribution of the excess energy Q cm determined with respect to the r\' meson 
production threshold for the pn reaction as obtained with sets of parameters given by the present 
calibration (solid histogram) and with the orientation of the spectator detector determined pre- 
viously from the analysis of the pn — > pnr\ reaction [42] (dashed histogram). 



As a cross check, the same data were analysed with the set of parameters describing 
the orientation of the spectator detector used previously in the analysis of the pn — > pnr\ 
reaction. The result is presented as a dashed histogram. 



6. Analysis of the experimental data 



In this chapter the determination of the four-momentum vectors of the outgoing nu- 
cleons will be presented, and the missing mass technique as well as the method of 
background subtraction will be described. 



6.1 Event selection 

The first step of the off-line analysis of the experimental data was the software selection 
of events with one track reconstructed in the drift chambers and a signal in the neutral 
particle detector. As a next step particles registered in drift chambers were identified. 

6.1.1 Proton identification 

All charged particles originating from the proton-neutron reaction are separated from 
the proton beam in the magnetic field of the dipole magnet due to the smaller momenta 
and some of them are registered in the drift chambers and scintillator detectors (SI, 
S3). For each event the time-of-flight of particle between the SI and S3 hodoscopes 
(At = ts3 — tsi) is measured. 




Figure 6.1: Distribution of the squared mass of charged particles originating from the pro- 
ton-deuteron reaction performed at a beam momentum of 3.35 GeV/c. 
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Knowing the relative distances between these detectors, the velocity of the particles 
is calculated. Based on the signals in the drift chambers Dl and D2 the trajectories 
of particles are reconstructed. Furthermore the particles momenta are reconstructed by 
tracking the trajectories back through the known magnetic field of the dipole to the 
centre of the interaction region (the overlap between the target and the beam). The pre- 
cision of the momentum reconstruction is equal to 6 MeV/c [54]. Having the velocity 
of the particle and its total momentum the squared invariant mass is determined. 
Charged ejectiles can well be identified as shown in figure 6.1, where the distribution 
of the squared mass of charged particles originating from the proton-deuteron reaction 
performed at a beam momentum of 3.35 GeV/c is presented. Three clear peaks are 
evidently visible and correspond to the squared mass of pions, protons and deuterons. 
For the further analysis events with invariant mass corresponding to the proton mass 
are selected (0.2 GeV 2 /c 4 < m 2 > 1.6 GeV 2 /c 4 ). 

After the particle identification, the time of the reaction (t r r eal ) at the target point is 
obtained from the known velocity, trajectory and the time measured by the SI detector. 

6.1.2 Neutron identification 

The determination of the four-momentum of the neutron is based on the time-of-flight 
measurement. The absolute value of the neutron momentum can be expressed as: 



The m n denotes the mass of the neutron, / is a distance from the centre of the interaction 
region to the central part of the module which fired as the first one. Thus the time-of- 
flight TOF N is taken as a difference between the time of the reaction f r eac (known from 
the backtracking protons to the interaction point) and the shortest time of the neutral 
particle detector. The individual momentum components are determined from the angle 
defined by the centre of the hit segment. The momentum component along the y-axis 
is assumed to be zero since the sensivity of the neutral particle detector is insufficient 
in the up-down direction. The granularity of the detector allows to determine the hori- 
zontal hit position with an accuracy of ±4.5 cm. 

Figure 6.2 presents the time-of-flight distribution between the target and the neutral 
particle detector as measured for neutral particles. The veto detector installed in front 
of the neutral particle detector discriminates signals originating from charged particles, 
however, it does not discern gamma quanta, which are detected in the neutral particle 
detector as well. Therefore, a clear signal originating from the gamma rays is seen at 
24.5 ns over a broad yield from neutrons. This histogram shows that a discrimination 
between signals originating from neutrons and gamma quanta can be done by a cut 
on the time-of-flight. From the Monte Carlo simulations of the pn — > pnr\' reaction 
the largest expected momentum value of the neutron is eqaul to 1.4 GeV/c which cor- 
responds to the time-of-flight value of 28 ns (see dashed histogram in fig. 6.2). This 




(6.1) 
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allows to select events for which the time-of-flight is higher than 26 ns as it is indicated 
by an arrow in fig. 6.2. This condition excludes background of gamma quanta. 
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Figure 6.2: Distributions of the time-of-flight between the target and the neutral particle detec- 
tor (T N ) for neutral particles. The experimental distribution depicted as a black solid histogram 
is compared with i) a simulated spectrum of the pn — ► pnr\' reaction (dashed histogram) and if) 
simulated exemplary background contributions (dotted histograms). 

Background reactions with neutrons in the exit channel and their contribution to the 
time-of-flight distribution have been investigated also. The simulated reactions were 
analysed in the same way as the experimental data. The distributions of the time-of- 
flight for neutrons originating from quasi-free pn — > pn pions reactions are shown in 
figure 6.2 with dotted histograms. The time-of-flight cut does not exclude neutrons 
from different quasi-free reactions, and this will constitute the main background which 
must be subtracted (see section 6.3.1). 

6.2 Determination of the excess energy 

Having proton and neutron identified, the consecutive software selection of the experi- 
mental data was applied. In addition to one proton registered in the drift chambers and 
a neutron in the neutral particle detector, at least a signal in one layer of the spectator 
detector was demanded. 

As it was disscused in chapter 3, to establish the exces energy Q cm with respect to 
the pn — > pnr\' reaction, the measurement of the proton spectator momentum, and the 
determination of the total energy y/s is indispensable. 



EXP : pn -> pn ij' 
SYM : pn -> pn r\' 
pn -> pn 7i° 

pn -> pn 2ji' 
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Figure 6.3: Simulated distribution of the Fermi momentum of nucleons inside the deuteron 
based on the Paris potential [60]. Arrows depict values of momenta for which the spectator 
protons are stopped in the first or in the second layer of the spectator detector. 



The spectator detector is designed such, that protons impinging perpendicular to the 
detector surface and having a kinetic energy of less than 6 MeV are stopped in the first 
layer (dE). Protons with kinetic energy in the range between 6 MeV and 9 MeV are 
stopped in the second layer (E). All protons with T spec > 9 MeV are passing through 
both layers of the spectator detector. As it is seen in figure 6.3 a large number of spec- 
tator protons (with maximum Fermi momentum of 130 MeV) will be stopped in the de- 
tector. In the experiment we registered the spectator protons with momenta larger than 
35 MeV/c since it is not possible to separate the energy losses for protons with lower 
momenta from the noise level. The spectator noise cut has been performed for each pad 
of the detector separately, using the spectra of energy loss triggered by a pulser with a 
frequency of 1 Hz additionally to other experimental triggers. Figure 6.4 (left) shows 
the energy losses in the I st layer of the spectator detector versus the 2 nd layer. Slow 
spectator protons are stopped in the first or second layer of the detector whereas fast 
particles (mainly charged pions) cross both detection layers. Signals from deuterons 
are not observed since deuterons cannot be emitted backward due to the kinematics. 
Having the deposited energy (T spec ) and the emission angles (Q S pec,§spec) we calculate 
the energy of the spectator proton (E spec = T spec + m p ) and its momentum (p sp ec)- Fig- 
ure 6.4 (right) presents the momentum distribution of protons considered as spectator 
(points) compared with simulations taking into account a Fermi motion of nucleons 
inside the deuteron (solid line). 
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Figure 6.4: Left: Energy losses in the first layer versus the second layer as measured with a 
deuteron target and a proton beam with a momentum of 3.35 GeV/c. Right: Momentum dis- 
tribution of the spectator proton as reconstructed in the experiment (points) in comparison with 
simulation calculations taking into account Fermi momentum distribution of nucleons inside 
the deuteron according to the Paris potential [60] (solid histogram). The simulated events were 
analysed in the same way as the experimental data. 

The total energy (\/s) available for the quasi-free proton-neutron reaction can be 
calculated for each event from the vectors of the spectator momentum and beam mo- 
mentum. 

s = (E b + E t f - (p b + p t ) 2 , 
Ef = m c i E speCl 

Pt = -Pspec (6.2) 

Eb and p b are the energy and momentum vector of the proton beam. E t and p t are 
the energy and momentum vector of the neutron target, calculated from the energy and 
the momentum of spectator proton (E spec , p spec )\ and m p denote the deuteron and 
proton masses, respectively. The excess energy with respect to the pn — > pnr\' process 
is equal to: 

Qcm = Vs~ (m p + m n + m^), (6.3) 

where m p , m n and m^i are the mass of proton, neutron and the T|' meson, respectively. 
The experimental distribution of the excess energy is shown in the left panel of fig- 
ure 6.5. As was previously discussed, with one fixed value of the beam momentum - 
due to the Fermi motion of nucleons inside the deuteron - one can scan a broad range 
of the excess energy, from up to 32 MeV for pbeam = 3.35 GeV/c. Negative values 
of the excess energy correspond to events, when the total energy in the centre-of-mass 
frame (y/s) was insufficient for the production of a T|' meson and only the production 
of multi-pion background events could occur. For excess energies larger than zero, ad- 
ditionally to the pion production also the Tj' meson is created. A simulated distribution 
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of the excess energy for quasi-free pn — > pnr{ reaction, calculated with respect to the 
pnr\' system is depicted in figure 6.5 (right). 
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Figure 6.5: Distributions of the excess energy Q cm for the quasi-free pn — > pnX reaction deter- 
mined with respect to the pnr\' threshold as obtained experimentally (left) and for Monte Carlo 
simulation (right). In the simulation all known detector features were taken into account (for 
details see text). 



A number of factors influence the registered yield of the pn — > pnr\' reaction as 
a function of the excess energy: i) the Fermi momentum distribution of nucleons in 
the deuteron target, ii) the relative settings of the spectator detector and the target, 
Hi) the absolute momentum of the proton beam and iv) the acceptance of the COSY-1 1 
detection setup. Due to the rapid growth of the total cross section close to the production 
threshold the accuracy of the excess energy determination is one of the most important 
parameters and it is mandatory to estimate correctly its uncertainty. The accuracy of 
the excess energy determination was estimated based on Monte Carlo simulations. For 
this purpose we have simulated Nq = 10 9 events of the pn — > pnr\' reaction taking 
into account the Fermi motion of the nucleons inside the deuteron, the size of the target 
(diameter ~ 9 mm), the spread of the beam momentum (opb = 1.5 Me V/c) and the 
horizontal (o(x) ~ 2 mm) and vertical (o(y) ~ 4 mm) beam size [101]. Simulated data 
have been analysed in the same way as the experimental data, taking into consideration 
the energy- and time-resolution of the detectors. 
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Figure 6.6: Difference between the generated (£2fm ) and reconstructed {ffcni) excess energy for 
the pn — > pnr\' reaction simulated for a proton beam momentum of 3.35 GeV/c and a deuteron 
target. 

Figure 6.6 presents the distribution of differences between the generated excess en- 
ergy (Qan) and the reconstructed excess energy from signals simulated in the detectors 
(Qan)- The distribution of A Q cm = Qan — 2^ was fitted by a Gaussian function 
resulting in an excess energy resolution of o(Q cm ) = 4.2 MeV. For the pn — > pnr\ 
reaction near threshold a standard deviation of the excess energy was derived to be 
G(Qcm) = 2.2 MeV [42], which is comparable with a (Q cm ) = 1.8 MeV obtained at 
similar conditions at the PROMICE/WASA setup [80]. 

6.3 Identification of the pn —> pnr\' reaction 

The Tj' mesons produced in the pn — > pnr\' reaction are identified using the missing 
mass technique. 

Knowing the four-momentum vectors of a proton and a neutron in the initial and fi- 
nal state, and employing the principle of momentum and energy conservation one can 
calculate the squared mass of the unmeasured particle or group of particles: 

m\ = E 1 x - p x 2 = (E b +E t —E p — E n ) 2 — (p b + Pt-p P ~Pn) 2 (6.4) 

where, 

Etnpb is the energy and momentum of proton beam, 
E t , p t is the energy and momentum of neutron target, 
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E p ,p p is the energy and momentum of outgoing proton, and 
E n ,p n is the energy and momentum of outgoing neutron. 

The method of the proton identification and of the neutron selection was described in 
sections 6.1 and 6.2. 

Before the missing mass was calculated, the collected data were grouped according to 
the excess energy. The available range of the excess energy Q cm above the r\' meson 
production threshold (from to 32 MeV) has been divided into four bins, each 8 MeV 
wide. The choice of the interval width is a compromise between the statistic and the 
resolution and reflects the accuracy (FWHM) of the determination of the excess energy. 
Due to the very low statistic and small acceptance, events corresponding to the excess 
energy interval between 24 and 32 MeV are not taken for further analysis. Moreover, 
in this range the excess energy distribution is very sensitive to changes of the position 
of the spectator detector, what can be seen from figure 6.7 showing the distribution 
of the excess energy Q cm determined for the pn — > pnr\' reaction as obtained with the 
position of the spectator detector fixed based on the data (solid histogram) and with the 
detector shifted from this position in the z-axis direction by 3 mm (dashed histogram). 
A significant reduction of number of counts in the interval between 24 and 32 MeV is 
clearly observed. 




Figure 6.7: The distribution of the excess energy Q cm determined for the pn — ► pnr\' reaction 
as obtained with the position of the spectator detector fixed based on the data (solid histogram) 
and with the detector shifted from this position in z-axis direction by 3 mm (dashed histogram). 
The shift of 3 mm is done for demonstration purpose. Note that in chapter 5.5.2 the accuracy of 
the position determination of the spectator detector was established to be ± 1 mm. 
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Figure 6.8: Experimental missing mass distri- 
butions of the pn — ► pnX reaction for ranges of 
the excess energy Q cm , from to 8 MeV (up- 
per panel left), from 8 MeV to 16 MeV (up- 
per panel right) and from 16 MeV to 24 MeV 
(lower panel). Arrows depict the nominal value 
of the Tj' mass. 



Figure 6.8 presents the experimental missing mass spectra of the pn — > pnX process 
calculated for the excess energy ranges Q cm = [0,8) MeV Q cm = [8, 16) MeV and 
Qcm = [16,24) MeV. For all missing mass distributions we observe an increase of 
counting rate towards the higher masses ending at the kinematical limit. The shape is 
mostly due to the detector acceptance. A signal from T|' meson production is expected 
on top of the multi-pion mass distribution at the position corresponding to the T|' mass - 
denoted by the arrow - however no enhancement around this region is seen. Therefore, 
to extract number of T|' events, the multi-pion background must be disentangled from 
the missing mass distributions. 

6.3.1 Method of background subtraction 

The background investigation in the missing mass spectrum of the quasi-free pn — > pnX 
reaction is a challenging task, especially close to threshold. This is even more difficult 
when the resolution of the mass determination is comparable with the excess energy. 
The method of separating contributions from the multi-pion production in the missing 
mass spectrum of the quasi-free pn — > pnX reaction is described in details in a dedi- 
cated article [102]. Therefore, here only the general principle will be pointed out. 




60 



Analysis of the experimental data 



One can disentangle the number of pn — > pnr\' events from single- and multi-pion back- 
ground without conducting model-dependent simulations by comparison of the missing 
mass distributions for the negative values of Q cm , when only pions can be created and 
for Q cm larger than zero. 




Figure 6.9: Schematic presentation of the background construction and subtraction, a) Dis- 
tribution of the excess energy with respect to the pn — > pnr\' reaction, b) and c) Form of the 
missing mass spectra as derived for negative (Qb) and positive (Qr\>) values of Q. d) The shape 
of the background constructed from the shape of the multi-pion mass distribution (shaded his- 
togram), shifted to the kinematical limit (dotted line) and normalised at low masses (dash-dotted 
line). The picture has been adapted from reference [102] 

Let us assume that the number of registered events is very large and that the shape 
of the reconstructed invariant mass of the multi-pion background is independent of 
the excess energy Q. The form of the background could then be determined from the 
missing mass spectrum dN n /dm from an infinitesimal range at any negative value of 
Q. The shape of the background can be expressed in terms of a normalised function 
B(m^i + Q — m) of the difference between the kinematical limit (mw + Q) and the 
given mass m. Assuming infinite statistics, one could divide the positive Q range into 
very narrow subranges within which the resultant missing mass spectrum would be a 
sum of a signal from the meson (r\') and the form B. This situation is schematically 
depicted in figure 6.9. In order to derive a signal of the meson from the missing mass 
spectrum for positive Q it would be sufficient to subtract from this spectrum a missing 
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mass distribution determined for negative Q values after the shift of the latter to the 
kinematical limit (dotted line) and normalisation at the very low mass values where no 
events from the Tj' meson are expected (dashed-dotted line) [102]. At real experimental 
conditions, due to the finite statistics, the missing mass spectrum is studied in finite 
intervals of Q, what alters the form of the multi-pion background. Therefore one has to 
take as narrow intervals of Q as reasonable, taking into account the experimental reso- 
lution in Q and assuming that within this range the multi-pion background is constant. 
An natural choice is to take the width of the Q subranges equal to the FWHM of the 
resolution of the Q determination. 
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Figure 6.10: Upper: Missing mass distribution determined for the pn — > pnim reaction at 
the excess energy range Q cm = (-30,-20) MeV (dotted line) and for the excess energy range 
Q cm = (-20,-10) MeV (solid line). The excess energy was calculated with respect to the pnr\' 
threshold. The dotted histogram was shifted to the kinematical limit of the solid histogram. 
Lower: Difference between the spectra of the upper panel normalised to the solid line. 



In order to gain more certainty that the measured shape of the multi-pion back- 
ground does not alter significantly with changes of the excess energy, we have per- 
formed simulations of the background originating from two pion creation in proton- 
neutron collisions. Figure 6.10 presents the missing mass distribution of the two pion 
system produced via the pn — > pnX reaction. Events were simulated at the excess en- 
ergy intervals of Q cm = (-30,-20) MeV and Q cm = (-20,-10) MeV with respect to the 
pnr\' system. From the comparison of the form of both spectra one realises that there 
is no noticeable difference in the shape at the missing mass range of the T|' mass. The 
systematic error caused by the subtracion of the shifted and normalised histogram is 
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shown in the lower panel of figure 6.10. It is demonstrated that the fractional system- 
atic error due to the multi-pion background subtraction is in order of one per cent of the 
background value within the range of 200 MeV in the vicinity of the kinematical limit. 

6.3.2 Background subtracted missing mass distribution 

As was already discussed, events with excess energies larger than zero were divided 
into intervals of AQ cm = 8 MeV. For each interval the missing mass was calculated 
separately. Next, for each distribution of the missing mass spectra, a corresponding 
background spectrum has been constructed from events with Q cm belonging to the 
( — 30, — 10) MeV range applying the method described in the previous section. As a 
cross check, to gain more confidence about the applied procedure we have constructed 
alternative background distrubutions taking events from a Q cm range from -20 MeV to 
MeV. The result of this investigation has shown a consistency of both background 
shapes in the order of 1%. 
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Figure 6.11: Experimental missing mass dis- 
tributions of the pn — ► pnX reaction (solid 
black line). The distributions were obtained 
for ranges of the excess energy Q cm , from 
to 8 MeV (upper panel left), from 8 MeV to 
16 MeV (upper panel right) and from 16 MeV 
to 24 MeV (lower panel). The corresponding 
background spectra (dashed red line) were con- 
structed as described above. 



Afterwards, to each distribution for Q cm > the corresponding background spec- 
trum has been normalised for mass values smaller than 0.25 GeV/c 2 where no events 
from the r\' are expected. In this region events correspond to the one pion production, 
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for which the total cross section remains nearly constant in the shown range of the 
excess energy. Figure 6.11 shows the experimental mass distribution for signal (solid 
black histogram) as determined for the excess energy ranges [0, 8) MeV, [8, 16) MeV, 
and [16,24) MeV. The dashed red line depicts the background spectrum, shifted to the 
kinematical limit and normalised accordingly. 
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Figure 6.12: Experimental missing mass distri- 
butions of the pn — ► pnr\' reaction for ranges of 
the excess energy Q cm , from to 8 MeV (up- 
per panel left), from 8 MeV to 16 MeV (up- 
per panel right) and from 16 MeV to 24 MeV 
(lower panel) determined after the background 
subtraction. Vertical bars indicate statistical er- 
rors. Arrows depict the nominal value of the T|' 
mass. 
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Figure 6.12 shows distributions of the missing mass for the pn — > pnr\' reaction 
as determined after subtraction of the background. Due to the low statistics and very 
low signal-to-background ratio the signal from T|' meson created in the proton-neutron 
collision is statistically insignificant. Therefore, in this thesis we can only estimate the 
upper limit for the Tj' meson production in the pn — > pnr\' reaction. 

Expected missing mass distributions from the Monte Carlo simulations are pre- 
sented in figure 6.13. In the simulation, based on the GEANT-3 package we took into 
account the momentum spread of the beam, beam and target dimensions, Fermi motion 
of nucleons inside the deuteron target, the geometry of the COSY-11 detector setup, 
multiple scattering of particles and other known physical effects were included in the 
GEANT-3 package [96]. The simulated data have been analysed in the same way as 
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the experimental one, taking into consideration position-, time- and energy-resolution 
of all detector components. 
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Figure 6.13: Simulated missing mass distribu- 
tion of the pn — ► pnr\' reaction for three ranges 
of the excess energy Q cm , from to 8 MeV (up- 
per left), from 8 MeV to 16 MeV (upper right) 
and from 16 MeV to 24 MeV lower panel. 
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In order to estimate the systematic error of the number of the pn — > pnX events, the 
change of the missing mass spectra has been studied by varying different parameters 
describing the experimental conditions (position of the detectors, time offsets, etc.) 
in the analysis and in the simulation. Change of the global time offset of the neutral 
particle detector by one standard deviation of its time resolution (of = 0.4 ns) resulted 
in systematic error of the missing mass distribution of 3%. As was disscused above, 
the systematic error due to the method used for the background subtraction is equal to 
1%. The uncertainty due to the model dependence calculation of the Fermi momentum 
distribution is equal to 2% [82] and was estimated as the difference between results 
determined using the Paris [60] and the CD-Bonn [103] potentials. A displacement 
of the spectator detector (fixed based on the experimental data) by 1 mm (standard 
deviation of its position accuracy) changes the number of events by 5%. 




7. Luminosity determination 



For the calculation of the absolute total cross section of the pn — > pnr\' reaction the 
knowledge of the integrated luminosity is mandatory. The determination of the lu- 
minosity is based on the registration of the quasi-free pp — > pp reaction, which was 
measured simultaneously with the r\' meson production in proton-neutron collisions. 
In this kind of scattering, the proton from the beam interacts with the proton bound 
inside the deuteron target. The outgoing protons were measured in coincidence with 
the SI and S4 scintillator detectors (see figure 4.2). In the S4 scintillator and the posi- 
tion sensitive silicon detector - referred to as the monitor detector - the recoil proton is 
registered. The forward proton is bent in the magnetic field of the dipole towards the 
drift chambers and SI scintillator detector. Figure 7.1 shows the measured correlation 
between the position in the SI scintillator and the position in the monitor detector. The 
lower band arises from the free proton-deuteron scattering, whereas the broad distribu- 
tion results from the quasi-free proton-proton scattering events. 
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Figure 7.1: The hit position of the recoil particles measured in the monitor detector versus the 
position in the SI detector of forward scattered particles. The pads number 15, 32 and 39 have 
lower counting rate due to a failure in their signal amplification. This imperfection was taken 
into account in the analysis. 

In the off-line analysis we have separated the quasi-free proton-proton from the 
proton-deuteron scattering applying a cut on the histogram as depicted by black line in 
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figure 7.1. For the further analysis only events above the black line were taken. The 
determination of the luminosity by means of the quasi-free proton-proton reaction mea- 
surement takes advantage of the availability of the precise cross sections determined by 
the EDDA group [104]. The normalisation uncertainty of the EDDA differential cross 
sections is relatively small (equal to circa 4% [104, 105]), whereas the systematic error 
arising from the determination of the corresponding differential cross section for the 
proton-deuteron scattering is in order of 10% [106]. 

The momentum of the fast proton is determined by tracking back the trajectory recon- 
structed from the signals in the drift chambers through the magnetic field to the target 
point. Figure 7.2 (left) shows the parallel versus the transversal component of the re- 
constucted momentum of the forward scattered proton. Events corresponding to the 
elastically scattered protons are seen near the kinematical ellipse, which is marked as 
a dashed line. The right panel of figure 7.2 presents an analogous spectrum obtained 
from the Monte Carlo simulations taking into account the acceptance of the detection 
system, the Fermi motion of the nucleons, and the variation of differential cross section 
for the elastic scattering as a function of the scattering angle and energy as it will be 
described later in this section. 
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Figure 7.2: Transversal versus parallel momentum component of the reconstructed forward 
proton momentum as obtained in the experiment (left) and in the simulation (right) 



In proton-proton experiments the elastic scattering pp — > ppis taken as a reference 
reaction due to the wealth of high-precision data available from the measurements of 
the EDDA collaboration [104]. The additional advantage of selecting this channel is 
the relatively high number of events, being clearly selectable by a coincident detection 
of both ejectiles. To determine the luminosity, the known differential cross sections for 
elastically scattered protons are compared to the number of scattered protons from the 
experimental data. 

In case of quasi-free elastic scattering we have to deal with the Fermi motion of the 
nucleons inside the deuteron. This motion implies that the value of the total energy in 
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the centre-of-mass system as well as the direction of the centre-of-mass velocity varies 
from event to event. To demonstrate the momentum spread caused by the Fermi motion 
let us consider the equivalent effective beam momentum as it is seen from the proton 
inside the deuteron. Figure 7.3 shows the distribution of the effective beam momentum 
for the quasi-free proton-proton scattering at the nominal value of the beam momen- 
tum equal to 3.35 GeV/c. The momentum changes significantly from 2.2 GeV/c up to 
4.5 GeV/c, therefore this effect cannot be neglected. At this point we have to emphasise 
that due to this effect, in a single subrange of the scattering angle in the laboratory sys- 
tem there are events originating from scattering at different values of the total energy 
y/s, as well as different scattering angles in the proton-proton centre-of-mass system. 
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Figure 7.3: Distribution of the effective beam momentum for the quasi-free proton-proton scat- 
tering at the nominal proton beam momentum of pbeam = 3.35 GeV/c as it is seen from the 
proton in the target. 

Therefore in order to calculate the integral luminosity we have to perform simula- 
tions taking into account these effects. Each simulated event has been associated with 
a weight corresponding to the differential cross section which is a function of the scat- 
tering angle and the total energy in the proton-proton centre-of-mass system \fs [107]. 
For this purpose we have used cross section values for the pp — > pp reaction computed 
by means of the SAID programme [108] because the accessible data base of the EDDA 
collaboration was insufficient. As it is seen in figure 7.3 the effective beam momentum 
which is seen from the nucleon inside the deuteron changes from 2.2 GeV/c up to 4.2 
GeV/c whereas EDDA measurements were performed in the beam momentum range 
from 0.712 GeV/c to 3.387 GeV/c. 
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In order to calculate the integrated luminosity, available range of the x-coordinate 
along the SI detector has been divided into four subranges. Further on, in order to 
separate the background originating from multi-particle reactions, the projection of the 
distance of the points from the kinematical ellipse was extracted for each subrange 
separately. The results for two subranges are shown in figure 7.4. Next, after the 
background subtraction the real number of scattered events AN eX p(Qiab) mt0 a given 
subrange of the SI detector has been obtained. 




distance to ellipse [a.u.] distance to ellipse [a.u.] 

Figure 7.4: Projection of the experimental event distribution from fig. 7.2 (left) along the ex- 
pected kinematical ellipse for two subranges of the SI detector. The dash lines depict the 
background as interpreated under the assumption that it is linear. The solid lines outline the 
background described by an exponential function. 




Figure 7.5: Projection of the simulated event distribution from fig. 7.2 (right) onto the expected 
kinematical ellipse for two subranges of the S 1 detector. 
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In order to determine the ANMc(®iab) number, we have simulated iVo = 10 8 quasi- 
free pp — ► pp events, which have been analysed using the same procedure as in case 
of the experimental data and which were weighted with the values of the differential 
cross sections expressed in units of mb. The results for analogous subranges as for the 
experimental distributions are shown in figure 7.5. 

The integrated luminosity was calculated as follows [107]: 

L= No ANexp(Qlab) = 

2%AN M c(Qiab) 
N()AN exp 

— 3£ (7.1) 

2n fAa(Q lab $ lab ) m ( *> PF, 9f, §F)f(pF, Qf, §F)dp F dcosd F d§ F d$*dcosQ* 

We have determined the integral luminosity for all subranges of the SI de- 
tector individually. The weighted average integrated luminosity equals to 
L = (4.77 ± 0.06) x 10 cm -2 . The estimated uncertainty is statistical only. 
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Figure 7.6: Differential cross sections as a function of the beam momentum for a few values 
of the scattering angle 6 cm in the centre-of-mass system. Black points stand for EDDA col- 
laboration data [104], lines denote SAID calculations [108]. Distribution of the effective beam 
momentum for quasi-free pp — > pp reaction calculated at the beam momentum of 3.35 GeV/c 
is also superimposed in the figure. 

In addition we have estimated the systematic uncertainty. One of the possible 
sources of the systematic error may be attached to the background subtraction. For 
the evaluation of this systematic uncertainty, we have subtracted the background in two 
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ways, namely in one case the background was described by an exponential function, 
in the latter case we assumed a linear background (compare the dash and solid lines in 
figure 7.4). Finally the systematic error due to the background subtraction estimated as 
the difference between the results determined with the two methods is not larger than 
±4%. 

The systematic error originating from the assumption of the potential model of the 
nucleon bound inside the deuteron is equal to about 2% [82]. Another source of 
the systematic error originates from the approximation of the differential cross sec- 
tion by the calculations using the SAID procedure [108]. The total cross section 
for the proton-proton scattering process has been measured up to beam momentum 
Pbeam ~ 3.4 GeV/c, whereas the effective beam momentum for the studied reaction 
ranges up to 4.2 GeV/c. Therefore for events with effective Pbeam > 3.4 GeV/c the dif- 
ferential cross sections were interpolated using the SAID programme. Figure 7.6 shows 
a comparison of the existing differencial cross section from the EDDA measurement 
and the SAID calculations. In the same figure the distribution of the effective beam 
momentum is also shown. In order to estimate this uncertainty, we have calculated the 
luminosity assuming a linear dependence between the J£ and the momentum for events 
with effective beam momentum larger than 3.0 GeV/c. Performed calculations show 
that the result obtained under this assumption differs from the result presented above 
by 5%. 
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Figure 7.7: Relation between the centre-of-mass scattering angle and the position in the SI 
detector for the quasi-free proton-proton scattering at pteam = 3.35 GeV/c as has been obtained 
from Monte Carlo simulations of the quasi-free pd — > ppn sp reaction. 



To reduce this systematic error as much as possible, events with such scattering 
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angles in the centre-of-mass system Q cm have been selected for which the differential 
cross sections calculated using the SAID programme are in agreement with the ^(p) 
distributions measured by the EDDA collaboration (see figure 7.6). Figure 7.7 shows 
the relation between the available range of the Q cm angles and the position in the SI 
detector for the quasi-free proton-proton scattering at p\,eam = 3.35 GeV/c as has been 
obtained in the Monte Carlo simulation. The rectangular area denotes the applied cut 
on the position in the SI detector, which restricts the Q cm angle to the range from 48 to 
69 degrees a good agrement between the SAID calculations and EDDA data is given. 
Summarizing, we have estimated the overall systematic error of the integrated luminos- 
ity to be 7%. 
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8. Excitation function 



In the following section we estimate the upper limit of the total cross section for the 
quasi-free pn — > pnr\' reaction as a function of the excess energy. Furthermore, an 
upper limit of the ratio R^t = c(pp^pprfj ^ e determined and compared to the 
analogous ratio calculated for the T| meson [12, 42]. 



8.1 Upper limit of the total cross section 

The total cross section for the pn — > pnr\' reaction was calculated using the formula: 

Tnr N^' 
o T0T {Q cm ) = T v , , (8.1) 

L 'int x c e// 

where N^' is the number of r\' mesons for a given excess energy Q cm , L int is the lumi- 
nosity integrated over the time of the measurement and E e ff is the detection efficiency 
including the geometrical acceptance of the detector system. 

The luminosity was established from the number of the quasi-free proton-proton scat- 
tering events (see chapter 7). The acceptance of the detector setup and efficiency was 
determined based on Monte Carlo studies. After subtracting the missing mass distribu- 
tions for the negative values of Q cm from the spectra for positive values of Q cm - due 
to the very low signal-to-background ratio - the signal from the T|' meson was found to 
be statistically insignificant. Therefore, we estimate the upper limit of the total cross 
section for the quasi-free pn — > pnr\' reaction. 

The number of the T|' mesons can be calculated as the difference between the number 
of events (N SIG ) in the peak for the signal (solid line in figure 6.11) and the number of 
events (N BACK ) of the background in the peak (dash line in figure 6.1 1): 

N r\' = N SIG _ N BACK. ^ = j N SIG + N BACK (8 2 ) 

The range for the integration has been chosen based on the simulation of the missing 
mass distributions (see figure 6.13). Assuming that in the experiment no T|' mesons are 
observed (N^ = 0) the oA^ 1 value was used to determine the upper limit of the total 
cross section for the pn — > pnr\' reaction at a 90% confidence level. The result is shown 
in figure 8.1 and in table 8.1. 
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upper limit of 


Qcm [MeV] 


o(pn — ► />nT|') 




at 90% CL [nb] 


[0,8) 


63 


[8,16) 


197 


[16,24) 


656 



Table 8.1: Upper limit of the total cross section for the pn — > /wr)' reaction as a function of the 
excess energy. The excess energy intervals correspond to the binning applied. 

The horizontal bars in figure 8.1 represents the intervals of the excess energy, for 
which the upper limit of the total cross section was calculated. 




10 20 30 40 50 

Q [MeV] 

Figure 8.1: Total cross sections for the pp — ► ppr\' reaction as a function of the excess energy 
(black open squares). Upper limit for the total cross section for the pn — ► pnr\' reaction as a 
function of the excess energy (red dots). 



The efficiency and the acceptance of the COSY-1 1 detector system was estimated 
based on Monte Carlo simulations. The error of the detection efficiency (aE e ff) is de- 
termined to be 3%. The main contribution of this systematic error comes from the un- 
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certainty in the efficiency in the neutron registration, as was discussed in section 5.3.4. 
The total systematic error of o T0T (Q cm ) was estimated as the quadratic sum of inde- 
pendent systematic errors originating from luminosity calculation (gL,- ;j/ ps 7%), de- 
termination of the efficiency (oE e ff ~ 3%) and estimation of the number of T|' mesons 
(oAf 1 ' ps 6%) and is equal to ps 10%. 

8.2 Ratio 

Figure 8.2 presents the upper limit of the ratio R^/ = p^Z^^y) °f me tota l cross 
section for the pn — > pnv( and pp — ► ppr]' reaction as a function of the excess energy 
(red crosses). Analogous ratios but for the T] meson are also shown as open squares. 
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Figure 8.2: Upper limit of the ratio (R^) of the total cross sections for the pn —* pnr\' and 
pp — > ^/?ri' reactions (red crosses) in comparison with the ratio (R^) determined for the r\ 
meson (open squares). The superimposed line indicates a result of the fit to the R^ data taking 
into account the final state interaction of nucleons [42]. 

The total cross section for the pp — > ppr\' reaction was measured in former exper- 
iments [13, 14, 15, 16, 17]. It reveals a strong excess energy dependence, especially 
very close to threshold and hence this dependence must be taken into account when 
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comparing to the results for the pn — > pnr\' reaction which were established for 8 MeV 
excess energy intervals. Therefore, for a given interval of excess energy, we have de- 
termined the mean value of the total cross section for pp — > ppr\' reaction using the 
parametrisation of Faldt and Wilkin fitted to the experimental data [44, 45] (for details 
see appendix B). 





upper limit of 




upper limit of 


Qan [MeV] 


o(pn — > pnx\') 


o(pp — > ppV[') [nb] 


Rry> 




at 90% CL [nb] 




at 90% CL 


[0,8) 


63 


19.6 


3.2 


[8,16) 


197 


68.7 


2.9 


[16,24) 


656 


122.7 


5.3 



Table 8.2: The columns of the table include from left to right: Excess energy range; Upper limit 
of the total cross section for the pn — > pnr\' reaction; Mean value of the pp — > ppr\' total cross 
section according to the parametrisation given in section B.2; Upper limit of the ratio R n > as a 
function of the excess energy. 

The determined pp — > ppr\' total cross section values together with the upper limit 
of the ratio R^> are presented in table 8.2. In case of the r\ meson, the ratio of the 
total cross sections for the reactions pn — > pnr\ and pp — > ppr\ was determined 
to be 7?r| Rs 6.5 at excess energies larger than 15 MeV [12], what suggests the 
dominance of the isovector meson exchange in the production mechanism, whereas 
close to threshold the ratio falls down to about 7?^ 3 [42]. This decrease may be 
explained by the different energy dependence of the proton-proton and proton-neutron 
final state interactions [44, 45]. For the Tj' meson the upper limit of the ratio for the 
excess energy range from MeV to 8 MeV is nearly equal to values of the ratio obtained 
for the T] meson, whereas for larger excess energies (from 8 MeV to 24 MeV) the ratio 
is notably lower. 



9. Comparison with model predictions 



The total cross section for the pn — > pnr\' reaction close to the threshold appeared to be 
too small to reveal a signal in the missing mass spectrum with the statistics achieved in 
the discussed experiment. 

The established upper limit does not allow to verify the scenario of the production via 
the possible gluonic component in the T)' wave function [9, 10]. However, the results 
also do not exclude this hypothesis and may be treated as a promissing finding in this 
context since the established upper limits for the ratio R^y are lower than the values of 
the ratio for the T| meson. 



Recently, X. Cao and X. G. Lee [7] performed an analysis of pN — > pNr\' and 
TtN — > Nr\' reactions within an effective Lagrangian approach, assuming that the 
Sn(1535) resonance is dominant in the r\' production. Numerical results show that 
the n exchange is most important in the pN — > pNr\' process and predict a large ratio 
of a(pn — > pnr\') to o(pp — > ppr\') equal to R-q' = 6.5 in the excess energy range from 
MeV to 160 MeV. Thus the authors assumed that the production of the T)' meson 
proceeds in the same way as the production of the T| meson, and as a consequence pre- 
dicted a R^t value equal to the experimentally established Rr, value. 
Yet, the experimentally determined upper limit of R^ being the result of this thesis is 
considerably lower than the R^/ = 6.5 predicted in reference [7]. Therefore, on this 
stage of knowledge we can exclude the hypothesis that the dominant production mech- 
anism of the T)' meson in nucleon-nucleon collisions is associated with the excitation 
and de-excitation of the Sn(1535) resonance. 

Another model has been elaborated by L. R Kaptari and B. Kampfer [6], considering the 
production of the T)' meson in the pp — > ppr\' and pn — > pnr\' reactions close to thresh- 
old by analysing the data within a covariant effective meson-nucleon theory by includ- 
ing meson and nucleon currents with resonances Sn(1650), Pn(1710) and Pi3(1720). 
In the framework of this model a reasonable agreement with data on pp — > ppr\' has 
been achieved with a contribution of the meson conversion currents. The total contri- 
bution of nucleonic and resonance currents was found to be small. 
The predictions for the pn — > pnr\' reaction, based on this model [6] are seen as a dotted 
line in figure 9.1, and are confronted with the upper limits of the total cross sections for 
the pn — > pnV[' reaction (left) and the upper limits of the ratio Rw (right). 
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Figure 9.1: (Left:) Upper limit of the total cross section for the pn — > pnr\' reaction in com- 
parison to calculations assuming that the r\' production mechanism is associated with me- 
son exchange [6]. (Right:) Upper limit of the ratio (/?«') of the total cross sections for the 
pn — ► pnrf and pp — > ppT[ f reactions (red crosses). The dotted line represents results of calcu- 
lations from [6]. 

The experimentally determined values of the total cross sections for the pn — > pri\\' 
reaction, as well as the results achieved for the ratio Rq' do not exclude the calculations 
of ref. [6]. 

For the pp — > ppT\' reaction, calculations of K. Nakayama and collaborators [41] 
yields the Sn(1650) resonance as the dominant contribution to the production current, 
whereas the Pn(1880) resonance, mesonic and nucleonic currents were found to be 
much smaller. The combined analysis of the pp — > ppr\' reaction and photoproduction 
reactions was crucial for these findings. The details of the excitation mechanism of 
the S 1 1 resonance, however, are not constrained by the currently existing data. As 
the authors pointed out [41], to learn more about the relevant excitation mechanism, 
observables other than the cross sections, such as the invariant mass distribution [48] 
are necessary. 



10. Summary and perspectives 



The main goal of this dissertation is the determination of the excitation function of the 
total cross section for the quasi-free pn — > pnr( reaction near the kinematical threshold. 
The motivation of the presented experiment was the comparison of the pp — > ppr\' and 
pn — > pnr\' total cross sections in order to learn about the production mechanism of the 
T|' meson in the channels of isospin 1=1 and 1 = and to investigate aspects of the 
gluonium component of the T|' meson. 

The experiment, which is described in this thesis has been performed at the cooler 
synchrotron COSY in the Research Centre Jiilich by means of the COSY-1 1 detector 
system. For the purpose of this experiment the standard COSY-1 1 detector setup was 
extended by the neutral particle and spectator detectors. The quasi-free pn — > pnr\' 
reaction has been induced by a proton beam with a momentum of 3.35 GeV/c in a 
deuteron target. All outgoing nucleons have been registered by the COSY-1 1 detec- 
tors, whereas for the r\' meson identification the missing mass technique was applied. 

The upper limit of the total cross section for the pn — > pnr\' reaction in the ex- 
cess energy range between and 24 MeV determined for the first time ever is pre- 
sented in this thesis. Having the total cross section for the pp — > ppr\' reaction in 
the same energy range, measured in previous experiments, the upper limit of the ratio 
Rri' = o(pn — > pnr\') I c(pp — > ppf}') has been estimated. The achieved result differ 
from the analogous ratio obtained for the T] meson. This allows to conclude that the 
production mechanism of these mesons in nucleon-nucleon collisions are different and 
the meson r\' is not dominantly created via the excitation of the Sn(1535) resonance as 
it is the case for the T] meson. 

It is also shown, that the theoretical predictions assuming the production mechanism to 
be associated with the fusion of virtual mesons or with glue excited in the interaction 
region of the colliding nucleons are not excluded by the experimental data. However, 
it is not possible to distinguish between these mechanisms based on the existing data. 
This conclusion calls for further theoretical as well as experimental investigations. 
For quantitative tests of these mechanisms an order of magnitude larger statistics and a 
larger energy range would be required. This can be reached with the WASA-at-COSY 
facility. 

In parallel to the investigation of the pn — > pnr\' reaction, presented in this thesis, 
we have extended the experimental studies to a pure isospin zero state of the interact- 
ing nucleons by the measurement of the quasi-free pn — > dr\' reaction. This experiment 
was conducted using the proton beam with a momentum of 3.365 GeV/c and a deuteron 
target. Assuming that the ratio of the cross sections for the pn — > dr\' and pp — > ppr\' 
reactions will be in the same order of magnitude as the ratio already established [12] for 
the pn — > dr\ and pp — > ppr\ reactions it is expected to identify about 1000 pn — > dr\' 
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events in the data sample [109]. 

Together with the results of the pp — > ppr\' and pn — > pnr\' reactions we would then 
complete the study of the r\' meson production cross section in nucleon-nucleon colli- 
sions. This will further reduce significantly the ambiguities of theoretical models and 
will lead to a better understanding of the production mechanism of the r|' meson in 
nucleon-nucleon collisions. The result will also be of importance for rate estimates for 
studying the r\' meson decays with the WASA-at-COSY facility [67, 110]. 



A. Precision of the neutron momentum 
determination 



The experimental precision of the missing mass determination of the pn — > pnX re- 
action strongly relies on the accurate measurement of the momentum of neutrons. The 
neutron detector is designed to deliver the information about the time at which the 
registered neutron or gamma quantum induced a hadronic or electromagnetic reaction. 
The time of the reaction combined with this information allows to calculate the time- 
of-flight (TOF^) of the neutron (or gamma) on the 7.36 m distance between the target 
and the neutron detector. Therefore the accuracy of the time-of-flight depends on a) 
the neutron detector time resolution and b) the time resolution of the S 1 counter. The 
neutral particle detector time resolution was determined to be 0.4 ns [91] which cor- 
responds to the time of flight in the scintillator material for a gamma quantum on a 
distance of a size of a single module. The neutron path (1) is defined as the distance be- 
tween the position of the centre of the target-beam overlap and the centre of the module 
which provided a signal as the first one. The accuracy of this measurement amounts 
approximately to half the size of the module which is 4.5 cm. In case of neutrons the 
absolute value of the momentum (p) can be expressed as: 



where m denotes the mass of the neutron, / stands for the distance between the target 
and the neutral particle detector and TOF N is the time-of-flight of the particle. Monte 
Carlo studies of the quasi-free pp — > pp reaction via the dp — > ppn spec process have 
been performed in order to establish the momentum resolution of the neutron detector. 
A spectator neutron has been denoted as n spec . This reaction was chosen since it could 
be identified in the data taken during previous COSY- 1 1 measurements conducted with 
the deuteron beam [111], and hence the simulation results could have been corrobo- 
rated by the experimental data [79]. The detection of both outgoing protons from the 
dp — > ppn sp reaction proceeds via the well established method [62] and the spectator 
neutron is measured by the neutral particle detector [79]. 

Figure A. 1 (left) presents the difference (AP = P gen — P rec ) between the generated 
neutron momentum (P gen ) and the reconstructed neutron momentum from signals sim- 
ulated in the detectors (P rec ). The value of (P rec ) was calculated taking into account 
the time resolution of the neutral particle detector (o = 0.4 ns) [91] as well as the time 
resolution of the SI counter (o = 0.25 ns) [77]. The time at which the proton hits the 
S 1 detector is needed to determine the time of the reaction in the target place. The dis- 
tribution of AP - for neutrons possessing a momentum of 1.6 GeV/c - was fitted by a 
Gaussian function resulting in a momentum resolution of o(P) = 0.14 GeV/c [79]. The 



/ 



1 



(A.1) 



p = m • 
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accuracy of the TOF N mesurement is approximately independent of the momentum of 
neutron, however, the relative resolution for determining the neutron momentum alters 
significantly as its momentum changes. This fractional momentum resolution - which 
is shown in figure A. 1 (right) as a function of the neutron momentum - is given by the 
equation: 

dp 

(A.2) 



P 



where o ? = y + a|j accounts for the time resolution of both the neutral particle 
and the SI detectors. The result presented in figure A.l was obtained assuming that 
o n = 0.4 ns, O51 = 0.25 ns, and / = 7.54 m. The value of / corresponds to the distance 
between the target and the centre of the neutron detector. 
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Figure A.l: Left: Difference between generated (P gen ) and reconstructed (P rec ) neutron mo- 
menta for the dp — > ppn spec reaction. Right: Fractional momentum resolution of the neutraf 
particfe detector as a function of neutron momentum. Arrows depict the kinetic energy and 
momentum range of neutrons originating from pd — > p sp pnX[' reaction at the beam momentum 
of 3.35 GeV/c. 



This fractional resolution depends on the neutron momentum and e.g. for the neu- 
trons produced at threshold for the pn — > pnr\ reaction (p = 0.76 GeV/c) it amounts to 
2% and at threshold of the pn — > pnr\' reaction (p = 1.06 GeV/c) it is equal to 3.2% 



B. Parameterisation of the pp — > ppr\ f 
total cross section 



In general the total cross section is defined as an integral of the probabilities to populate 
a given phase space interval over the whole kinematically available range normalised 
to the flux factor F according to formula: 



dV ps \M\ 



(B.l) 



where \M\ 2 denotes the production amplitude and V ps is the phase space volume. Fig- 
ure B.l shows the total cross section for the pp — > ppr\' reaction as a function of the 
excess energy. It is clearly seen that near thereshold the Tj' production strongly depends 
on the energy. 
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Figure B.l: Total cross sections for the pp — > ppr\' reaction as a function of the excess energy 
derived from the COSY-11, DISTO and SPESIII measurements [13, 14, 15, 16, 17]. The solid 
line depicts a parameterisation of the total cross section using formula B.2 



The solid line depicts a parameterisation of the total cross section utilizing a fac- 
torisation of the production matrix element \M\ into the short range primary amplitude 
and the initial and final state interaction. The proton-proton FSI has been taken into 
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account according to the Faldt and Wilkin model [44, 45] using the following formula: 

, . C 2 A? ( l + ^f)- 2 <B.2, 

y / X(5,mf,m^) 

where Ci, C2 and £ have to be determined from the data, and the triangle fuction 
X(x,y,z) [112] is defined as: 

X(x,y,z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz. (B.3) 

The function B.2 was fitted to the data and as a result values of the free parameters C2 
and £ have been obtained [113]: 

£ = (0.62 ± 0.13) MeV 
C 2 = (84 ± 14) mb 

Applying the parameterisation from equation B.2 with known values of C2 and £ one 
can compute the total cross section for the T|' meson production at any value of the 
excess energy close to threshold. 
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